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Based on the low energy effective Hamiltonian with generalized factorization, 
we calculate the new physics contributions to the branching ratios and CP-violating 
^3 ' asymmetries of the two-body charmless hadronic decays B — > PP, PV from the 

"fH , new strong and electroweak penguin diagrams in the Topcolor-assisted Technicolor 

CL|| (TC2) model. The top-pion penguins dominate the new physics corrections, and 

Q^' both new gluonic and electroweak penguins contribute effectively to most decay 

D . modes. For tree-dominated decay modes B — > vrvr, pir, etc, the new physics cor- 

. . i rections are less than 10%. For decays B — > K^*^Tr, K^*^r], ix^r}^ 

K K, etc, the new physics enhancements can be rather large ( from —70% to 
^ ■ ~ 200% ) and are insensitive to the variations of N^^f , A;^, rj and m^- within the 

' reasonable ranges. For decays B^ — > (p-K, 4>'i]^'\ K*K^ and p^K^, SB is strongly 

A^^'^'^— dependent: varying from —90% to ~ 1680% in the range of N^-f^l' = 2 — oo. 
The new physics corrections to the CP-violating asymmetries Acp vary greatly for 
different B decay channels. For five measured CP asymmetries of i? — > Kt:, Kr]' , ujtt 
decays, SAcp is only about 20% and will be masked by large theoretical uncertain- 
ties. The new physics enhancements to interesting B Kr]' decays are significant 
in size (~ 50%), insensitive to the variations of input parameters and hence lead 
to a plausible interpretation for the unexpectedly large B Kr]' decay rates. The 
TC2 model predictions for branching ratios and CP-violating asymmteries of all 
fifty seven B PP, PV decay modes are consistent with the available data within 
one or two standard deviations. 



PACS numbers: 13.25.Hw, 12.15.Ji, 12.38.Bx, 12.60.Nz 
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1 Introduction 



The main goals of B experiments undertaken by CLEO, BarBar, Belle and other collab- 
orations are to explore the physics of CP violation, to test the standard model (SM) at 
an unexpected level of precision, and to make an exhaustive search for possible effects of 
physics beyond the SM ^. Precision measurements of B meson system can provide 
an insight into very high energy scales via the indirect loop effects of new physics(NP). 
The B system therefore offers a complementary probe to the search for new physics at 
the Tevatron, LHC and NLC, and in some cases may yield constraint which surpass those 
from direct searches or rule out some kinds of NP models [|l|. 

In B experiments, new physics beyond the standard model may manifest itself, for 
example, in the following waysfl], |^: 

• Decays which are expected to be rare in the standard model are found to have large 
branching ratios; 

• CP-violating asymmetries which are expected to vanish or be very small in the 
SM are found to be significantly large or with a very different pattern with what 
predicted in the SM; 

• Mixing in B decays is found to differ significantly from SM predictions; 

These potential deviations may originate from the virtual effects of new physics through 
box and/or penguin diagrams in various new physics models @, |^, ||, |^, H, P]. 

Due to the anticipated importance of two-body charmless hadronic decays B ^ hih2 
( where hi and h2 are the light pseudo-scalar (P) and/or vector(V) mesons ) in under- 
standing the phenomenon of CP violation, great effort have been made by many authors 
[0, [n], m, |l3l It is well known that the low energy effective Hamiltonian is the basic 



tool to calculate the branching ratios and Acp of B meson decays. The short-distance 
QCD corrected Lagrangian at NLO level is available now, but we do not know how to 
calculate hadronic matrix element from first principles. One conventionally resort to the 



factorization approximation |]T5|. However, we also know that non-factorizable contribu- 
tion really exists and can not be neglected numerically for most hadronic B decay channels. 
To remedy factorization hypothesis, some authors |jT6|, |1^ introduced a phenomeno- 
logical parameter N^ff (i.e. the effective number of color) to model the non-factorizable 
contribution to hadronic matrix element, which is commonly called generalized factoriza- 
tion. 



On the other hand, as pointed by Buras and Silverstrini [T^], such generalization 
suffers from the problems of gauge and infrared dependence since the constant matrix fy 
appeared in the expressions of C^^^ depends on both the gauge chosen and the external 
momenta. Very recently, Cheng et al. studied and resolved above controversies on 
the gauge dependence and infrared singularity of C^^^ by using the perturbative QCD 



factorization theorem. Based on this progress, Chen et al. |T4] calculated the charmless 



hadronic two-body decays of and B^ mesons within the framework of generalized 
factorization, in which the effective Wilson coefficients C^^^ are gauge invariant, infrared 
safe, and renormalization-scale and -scheme independent. 
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On the experimental side, the observation of thirteen B PP, PV decays by CLEO, 



BaBar and Belle collaborations |19, 20, 21, 23, 23, Ej, 25] signaled the beginning of the 



golden age of B physics. For B — > Kn, tttt decays, the data are well accounted for in the 
effective Hamiltonian p7|, p8| with the generalized factorization approach |jl5|, |12|, [l^. For 
B — Kt]' decays, however, the unexpectedly large decay rate B{B Kr]') = (SO^g^ ± 
7) X 10^^ ||2^ still has no completely satisfactory explanation and has aroused considerable 
controversy |2^ . 

In this paper, we will present our systematic calculation of branching ratios and 
CP-violating asymmetries for two-body charmless hadronic decays B PP, PV (with 
charged Bu, neutral Bd mesons ) in the framework of Topcolor- assisted technicolor (TC2) 
model by employing the effective Hamiltonian with the generalized factorization. 
Since the scale of new strong interactions is expected around 1 TeV, the tree-level new 
physics contributions are strongly suppressed and will be neglected. We therefore will 
focus on the loop effects of new physics on two-body charmless hadronic B meson decays. 
We will evaluate analytically all new strong and electroweak penguin diagrams induced by 
exchanges of charged top-pions tt^ and technipions vrf and vr^ in the quark level processes 
b sV* with V = 'J, gluon, Z, and then combine the new physics contributions with 
their SM counterparts, find the effective Wilson coefficients and finally calculate the new 
physics contributions to the branching ratios and CP-violating asymmetries for all fifty 
seven decay modes under consideration. We will concentrate on the new physics effects 
on charmless B PP, PV decays and compare the theoretical predictions in TC2 model 
with the SM predictions as well as the experimental measurements. For the phenomeno- 
logically interesting B Krj' decays, we found that the new physics enhancements are 
significant in size, ~ 50%, insensitive to the variations of input parameters and hence lead 
to a plausible interpretation for the large B Krj' decay rates. 

This paper is organized as follows. In Sec. 2, we describe the basic structures of the 
TC2 model and examine the allowed parameter space of the TC2 model from currently 
available data. In Sec. 3, we give a brief review about the effective Hamiltonian, and then 
evaluate analytically the new penguin diagrams and find the effective Wilson coefficients 
C^^^ and effective numbers with the inclusion of new physics contributions. In Sec. 4 
and 5, we calculate and show the numerical results of branching ratios and CP- violating 
asymmetries for all fifty seven B PP, PV decay modes, respectively. We concentrate 
on modes with well-measured branching ratio and sizable yields. The conclusions and 
discussions are included in the final section. 



2 TC2 model and experimental constraint 

Apart from some differences in group structure and/or particle contents, all TC2 models 
ppj , have the following common features: (a) strong Topcolor interactions, broken 
near 1 TeV, induce a large top condensate and all but a few GeV of the top quark mass, 
but contribute little to electroweak symmetry breaking; (b) Technicolor [32| interactions 



are responsible for electroweak symmetry breaking, and Extended Technicolor (ETC) |Q 
interactions generate the hard masses of all quarks and leptons, except that of the top 
quarks; (c) there exist top-pions vr^ and vr*^ with a decay constant Fj^ ^ 50 GeV. In this 
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paper we will chose the well-motivated and most frequently studied TC2 model proposed 
by Hill |30[ as the typical TC2 model to calculate the contributions to the charmless 
hadronic B decays in question from the relatively light unit-charged pseudo-scalars. It is 
straightforward to extend the studies in this paper to other TC2 models. 

In the TC2 modelBO], after integrating out the heavy coloron and Z\ the effective 



four-fermion interactions have the form [^| 



^1 

K - ^ ) i^ibRbRipL 



(1) 



where k = {g'^/iir) cot"^ 9 and Ki = {gf / in) cot"^ 9' , and My is the mass of coloron V"" 
and Z' . The effective interactions of (|l|) can be written in terms of two auxiliary scalar 
doublets 01 and 02- Their couplings to quarks are given by 



Ceff = Xli'L<PltR + A2^L02&/ 



(2) 



where Xf = 4tt{k + 2ki/27) and A2 = 47r(K — ni/27). At energies below the Topcolor 
scale A ~ 1 TeV the auxiliary fields acquire kinetic terms, becoming physical degrees of 
freedom. The properly renormalized 0i and 02 doublets take the form 



TT 



H+ 



(3) 



where tt^ and 7f° are the top-pions, and are the b-pions, ht is the top-Higgs, and 
Fjf ^ 50GeV is the top-pion decay constant. 

From eq.(^, the couplings of top-pions to t- and b-quark can be written as pO| : 



ml 



i ttn" + i tubLii^ + i —tibRTi^ + h.c. 



(4) 



where m* = (1 — e)mj and ml ^ IGeV denote the masses of top and bottom quarks 
generated by topcolor interactions. 

For the mass of top-pions, the current 1 — a lower mass bound from the Tevatron data 
is mjf > 150GeV^^, while the theoretical expectation is fa (150 — SOOGeV ) ||30| . For 
the mass of b-pions, the current theoretical estimation is m^^o ~ m^o ~ (100 — 350)GeV 

For the technipions vrf and vr^, the theoretical estimations 
200GeVp^, p8[. The effective Yukawa couplings of ordinary 



and mfj 



2mf §6 



are m,ri > 50GeV and m. 



technipions vrf and tt^ to fermion pairs, as well as the gauge couplings of unit-charged 
scalars to gauge bosons 7, and gluon are basically model-independent, can be found 



in refs. jS^, ^ 



At low energy, potentially large flavor- changing neutral currents (FCNC) arise when 
the quark fields are rotated from their weak eigenbasis to their mass eigenbasis, realized by 
the matrices Ul,r for the up-type quarks, and by Dl r for the down- type quarks. When 
we make the replacements, for example. 



bL 
bR 



Dl^dL + Dl'SL + D'jPbL, 



(5) 
(6) 
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the FCNC interactions will be induced. In TC2 model, the corresponding flavor changing 
effective Yukawa couplings are 



(7) 



For the mixing matrices in the TC2 model, authors usually use the "square-root 
ansatz": to take the square root of the standard model CKM matrix {Vckm = U^Dl) as 
an indication of the size of realistic mixings. It should be denoted that the square root 
ansatz must be modifled because of the strong constraint from the data of — mixing 
35| , In TC2 model, the neutral scalars and can induce a contribution to 



the B^ ~ = d, s) mass difference p5|] 



AM 



mr 



12 F|m|„ 



(8) 



where Mb^ is the mass of Bg meson, Fb^ is the i?q-meson decay constant, Bb^ is the 



renormalization group invariant parameter, and 6bq 



For Bd meson, using the 
= 50GeV, .[b^^Fb, = 



data of AMb^ = (3.05 ± 0.12) x IQ-^^MeV^ and setting F^ 
200Me\^, one has the bound 6bd < 0.76 x 10"'' for nifjo < QQQGeV . This is an important 
and strong bound on the product of mixing elements D^^^.- pointed in |^|, if one 
naively uses the square-root ansatz for both Dl and Dji, this bound is violated by about 2 
orders of magnitudes. The constraint on both and from the data of 6 ^ 57 decay is 
weaker than that from the B^—B^ mixings |]3^. By taking into account above experimental 
constraints, we naturally set that = for i ^ j. Under this assumption, only the 



charged technipions vr^^ , vTg and the charged top-pions vr contribute to the inclusive 
charmless decays h sqq, dqq with q G {m, d, s} through the strong and electroweak 
penguin diagrams. 

In the numerical calculations, we will use the "square-root ansatz" for D^^ and D^^, i.e, 
setting D^^ = Vtdf^ and Dff = Vts/2, respectively. We also flx the following parameters 
of the TC2 model in the numerical calculation [} 



lOOGeV, m^g 



200GeK F^ = SOGeK F^ = 120GeK e = 0.05, 



(9) 



where F^^ and Fj^ are the decay constants for technipions and top-pions, respectively. For 
m^f, we consider the range of m^- = 200 ± 100 GeV to check the mass dependence of 
branching ratios and CP-violating asymmetries of charmless B decays. 



Effective Hamiltonian and Wilson coefficients 



We here present the well-known effective Hamiltonian for the two-body charmless decays 
B —>■ hih2. For more details about the effective Hamiltonian with generalized factorization 



for B decays one can see for example refs.|12, 14, 27, 28 



^From explicit numerical calculations in next section, we know that the new physics contributions from 
technipions irf and tt^ are much smaller than those from top-pion tt* within the reasonable parameter 
space. We therefore fix — lOOGeV and mTr^ — 200GeV for the sake of simplicity. 
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3.1 Operators and Wilson coefficients in SM 

The standard theoretical frame to calculate the inclusive three-body decays b sqq ^ is 
based on the effective Hamiltonian p8| , |12|, 



10 

i=3 



CgQg 



{Saqp)v-A{q(5hc 



V-A, 



Here the operator basis reads: 

Qi = {sq)v-A{qb)v~A, Q2 
with q = u and q = c, and 

Q3 = {sb)v-Aj2((l'l')v-A, Q4 = {Sabi3)v-Aj2(l'p(la)v-A 
Q5 = 



Qr 
Q9 



q' q' 

{sb)v-Aj2i(l'(l')v+A, Qe = {Sabf3)v~Aj2(l' fi^'a)v+A, 
q' q' 

3 3 

2(^b)v~AY.^<l'(^'l')v+A, Qs = ^isMv-Aj2^<l'(l'f3(l'a)v+A, 
q' q' 

3 3 

-{sb)v-AY.^l'(l'l')v-A, QlO = -iSabf})v-AY.^<l'('i' I3(la)v-A, 
<?' Q' 

lL.m,-s^a^^{l + ^,)T:pbpG%, 



(10) 

(11) 

(12) 
(13) 

(14) 

(15) 
(16) 



where a and /3 are the SUiVj color indices, T^^ 



, 8) are the Gell-Mann matri- 



ces. The sum over q' runs over the quark fields that are active at the scale jj, = 0(1711,), 
i.e., q' G {u,d, s,c,b}. The operator Qi and Q2 are current-current operators, Q3 — Qq 
are QCD penguin operators induced by gluonic penguin diagrams, and the operators 
Qr " Qio are generated by electroweak penguins and box diagrams. The overall factor 
2/3 is introduced for convenience, and the charge e^/ is the charge of the quark q' with 
q' = u, d, s, c, b. The operator Qg is the chromo-magnetic dipole operator generated from 
the magnetic gluon penguin. Following ref.|T2|, we also neglect the effects of the electro- 
magnetic penguin operator Qj-y, and do not consider the effect of the weak annihilation 
and exchange diagrams. 

Within the SM and at scale Mw, the Wilson coefficients Ci{Mw), ■ ■ ■ , Cio{Mw) and 
Cg{M\Y) have been given for example in [ P7| , E 
regularization (NDR) scheme 



They read in the naive dimensional 



Ci{Mw) 
C2{Mw) 
Cs{Mw) 



11 as{Mw) 35 cte 



6 An 
n as{Mw) 
2 



18 471 



An 
as{Mw] 



247r 



Eo{xt) - - 



+ 



a. 



6n sin^ 9]y 

^For b — )■ dqq decays, one simply make the replacement s — > d. 



6 



Cg{Mw 

Cio{Mw 
Cg{Mw 



w 



Stt 
247r 



Eo{xt) - - 



in 



Eo{xt) - - 



0, 



Do{xt) - - 



Git 
0, 

E'oixt) 



4Co{xt) + Doixt) - - + 
y 



sin^ 9w 



-{lOBoixt) - ACoixt)) 



(17) 
fl8) 



rrif/M^^, the functions -80(3^)5 Co{x), Dq{x), Eo{x) and E'q{x) are the famihar 



where Xt 
Inami-Lim functions [E3 



Bo{x) 
Co{x) 
Doix) 
Eo{x) 



E'oix) 



1 
4 

X 



X 



+ 



x\nx 

1 — X {x — ly 

X — Q 3x + 2 
X — 1 (x — 1)^ 

-19x^ + 25x2 x^{5x^ -2x -6) 



\nx 



--lnx + 

9 36(a;-l)3 



+ 



18a; — llx^ — x^ 
12fl -x)3 



18(x - 1)4 
4 - 16a; + 9x2 

In X , 



Inx . 



2x + 5x 



4(1 -x) 



+ 



6(1 -x)4 
3x2 



2(1 



X 



log[x] 



(19) 
(20) 
(21) 
(22) 
(23) 



Here function Bo{x) results from the evaluation of the box diagrams with leaving lepton 
pair uu or /+/^|2^, function Co(x) from the Z'^-penguin, function Do{x) and Eo{x) from 
the photon penguin and the gluon penguin diagram respectively, and finally function 
Eq{x) arise from the magnetic gluon penguin. 

By using QCD renormalization group equations p7|, p^ , it is straightforward to run 
Wilson coefficients Ci{Mw) from the scale /i = O(Miy) down to the lower scale /i = 0(mb). 
Working consistently to the next-to-leading order ( NLO ) precision, the Wilson coeffi- 
cients Ci for i = 1, . . . , 10 are needed in NLO precision, while it is sufficient to use the 
leading logarithmic value for Cg. At NLO level, the Wilson coefficients are usually renor- 
malization scheme(RS) dependent. In the NDR scheme, by using the input parameters 
as given in Appendix A and setting /x = 2.5 GeV, we find: 



C5 
C9 



1.1245, C2 = -0.2662, C3 = 0.0186, C4 



-0.0458, 



0.0113, Cg = -0.0587, C7 = -5.5 x 10-^ Cg = 6.8 x 10" 
-0.0095, Cio = 0.0026, Cf^ = -0.1527. 



(24) 
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Here, 



Co + C*5- These NLO Wilson coefficients are renormalization scale and 



scheme dependent, but such dependence will be cancelled by the corresponding depen- 
dence in the matrix elements of the operators in Tie//? as shown explicitly in 



3.2 New strong and electroweak penguins in TC2 model 

For the charmless hadronic decays of B meson under consideration, the new physics will 
manifest itself by modifying the corresponding Inami-Lim functions Co{x) , Dq{x) , Eq{x) 
and E'q{x) which determine the coefficients Cs^Mw) , ■ ■ ■ , Ciq^Mw) and Cg^Mw), as il- 
lustrated in Eqs. (p!7| , p!8D . These modifications, in turn, will change for example the stan- 
dard model predictions for the branching ratios and CP-violating asymmetries for decays 
B hih2. 

The new strong and electroweak penguin diagrams can be obtained from the cor- 
responding penguin diagrams in the SM by replacing the internal lines with the 
unit-charged scalar {Tif,7if and vr^ ) lines, as shown in Fig.|l]. In the analytical calcula- 
tions of those penguin diagrams, we will use dimensional regularization to regulate all the 
ultraviolet divergences in the virtual loop corrections and adopt the MS renormalization 
scheme. It is easy to show that all ultraviolet divergences are canceled for each kind of 
charged scalars, respectively. 

we calculate analytically the new Z^- 



Following the same procedure of refs.[41, 43 



penguin diagrams induced by the exchanges of charged scalars vrf , vr^ and tt^, we find the 
new Co function which describe the NP contributions to the Wilson coefficients through 
the new Z°-penguin diagrams. 



r<TC'2 
^0 



1 



V2GfM^ 



m. 



mz 



4F. 



TTg 



with 



To(x) 



X 



X 



3Ff 



x] 



XI 



\og[x] 



(25) 



(26) 



where yt = iril^ /m\ with = (1 — e)mt, Zt = {^fntY / fn^-^i^t = (c'^t)^/''^^8- 

By evaluating the new 7-penguin diagrams induced by the exchanges of three kinds 
of charged pseudo-scalars (vf^^, vrj'', vr^), we find that. 



JJTC2 



-Mvt) + 



with 



Fn(x) 



47 - 79x + 38x^ 3 - 6x^ + Ax^ 



[Fo(^t) + 8Fo(6)] , 



log[x]. 



(27) 



(28) 



108(1 -x)3 18(1 -x)4 

By evaluating the new gluon-penguin diagrams induced by the exchanges of three 
kinds of charged pseudo-scalars (tt^*", nf,7cf) as being done in @, ^, we find that, 

1 ... 1 



E 



\4V2GfF? 

,TC2 _ j 1 



3V2GfF^ 
1 



6V2GfF^ 



[/o(;2i) + 8/o(6)+9iVo(6)]|, (29) 
[Ko{zt) + 8Ko{^t)+9Lom], (30) 
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with 



= 36(1 -x)3 -6(r^^°gt^]' ^''^ 
, , 5 - 19a; + 20^^ a;^ - 2x3 ^ ^ 

'^°(^' = - 6(1 +(r3i;)7'°^M- (32) 

^°<"' = - '^6(1-".^ -flTfj7'°sN. (33) 

- "36(1 -t)^ ^ 6(1^ (3^' 

Using the input parameters as given in Appendix A and Eq.(|^), and assuming m^f = 
200GeV, we find numerically that 

{Co, /^o, ^0, E'^Y^'X^M^ = {1.27, 0.27, 0.66, -1.58} (35) 

if only the new contributions from top-pion penguins are included, while 

{Co, /^o, ^0, E'^Y'^X^My, = {0.0002, 0.03, 0.04, -0.14} (36) 

if only the new contributions from technipion penguins are included. It is evident that 
it is the charged top-pion tt^ that strongly dominate the NP contributions, while the 
technipions play a minor rule only. We therefore fix the masses of and tt^ in the 
following numerical calculations. 

Using the input parameters as given in Appendix A and Eq. (H) and assuming m^f = 200 
GeV, we find that 

{Co,Do,E^,E'^f\^=M^ = {0.81,-0.48,0.27,0.19}, (37) 
{Co,Do,^o,^or''\=M^ = {1.27,0.30,0.71,-1.72}. (38) 

It is easy to see that the new physics parts of the functions under study are comparable 
in size with their SM counterparts. The SM predictions, consequently, can be changed 
significantly through interference. For Co and £^0 functions, they will interfere construc- 
tively. For Dq and Eq functions, on contrary, they will interfere destructively. One also 
should note that the magnitude of E''^^'^ is much larger than its SM counterpart, and 
hence E''^^^ will dominate in the interference. We will combine the two pats of the 
corresponding functions to define the functions as follows, 

Co(Mh/) = C,{Mwf^' + C,{Mwf''\ 

Dq{Mw) = D^iMw^^ + Dq{MwT''\ 

Eo{Mw) = E^{Mwf^' + Eo{Mwf''\ 

E'q{Mw) = E'^iM^f"' + E'q{MwY^\ (39) 

where the functions Dq{MwY^^ , Eo{Mw)^^^, Co{Mwf^^ and E'^iMwf^ have been 
given in Eqs.(|0,|l|, |2||1), respectively. While the functions Co(MvK)^^^ Doi^Mw f^'^, 
Eo{Mwf^^ and E'^^MwY^^ have also been defined in Eqs.(g|,g3,^,|0|), respectively. 
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Since the heavy charged pseudo-scalars appeared in TC2 model have been integrated 
out at the scale Mw, the QCD running of the Wilson coefficients Ci^Mw) down to the 
scale n = 0{miy) after including the NP contributions will be the same as in the SM. In 
the NDR scheme, by using the input parameters as given in Appendix A and Eq.(^), and 
setting m^j- = 200 GeV and /i = 2.5 GeV, we find that: 

Ci = 1.1245, C2 = -0.2662, C3 = 0.0195, = -0.0441, 

C5 = 0.0111, Cg = -0.0535, C7 = 0.0026, Cg = 0.0018, 

Cg = -0.0175, Cio = 0.0049, C^-^^ = 0.3735, (40) 

where Cg^^ = Cg + C^. By Comparing the Wilson coefficients in Eq. (^0|) with those given 
in Eq.(^3|), we find that Ci^2 remain unchanged, 6*3,4,5,6 changed moderately, 6*7,8,9,10 and 
Cg-^^ changed significantly because of the inclusion of new physics contributions. 



3.3 Effective Wilson coefficients 



Using the generalized factorization approach for nonleptonic B meson decays, the renor- 



malization scale- and scheme-independent effective Wilson coefficients C, 



10) 



have been obtained in 16 



|l^, I13] by adding to C*j(/i) the contributions from vertex- type 
quark matrix elements, denoted by anomalous dimensinal matrix 7y and constant matrix 
fv as given for example in [jl2l. Very recently, Cheng et al. [|18| studied and resolved the 
so-called gauge and infrared problems [0 of generalized factorization approach. They 
found that the gauge invariance is maintained under radiative corrections by working in 
the physical on-mass-shell scheme, while the infrared divergence in radiative corrections 
should be isolated using the dimensional regularization and the resultant infrared poles 



are absorbed into the universal meson wave functions ^ ^ 

In the NDR scheme and for S'f/(3)c, the effective Wilson coefficients C^^^ can be 



written as 12, 14 



(jeff 



fy + 7y log 



[fy+Jy log 



fy + 7y log 



fy + 'Jy log 



fy + Jy log 
fy + 'Jy log 



fy + 7y log 



/i J 
/i ; 

/i / 



6i 



2j 



3i 



4i 



5i 



6 47r 

2 47r 

6 An 
1 a„ 



(Ct + Cp + Cg) , 
(Ct + Cp + Cg) , 

(a + Cp + Cg) , 



Cj + ^^{Ct + C, + Cg) 



a. 



7i 



' C 



10 





Go 




- r + 


•-"10 





+ 7y log 



+ 7y log 



/i J 
/i J 



C. 



] 1 



8j 



Up 



9i 



-a 



+ 7y log 



(41) 



lOi 



where the matrices ry and 7y contain the process-independent contributions from the 
vertex diagrams. Like ref. |T^, we here include vertex corrections to Cj — Cio0. The 
anomalous dimension matrix 7^ has been given explicitly, for example, in Eq.(2.17) of 
T^ . Note that the correct value of the element {ffqDR)^^ and {fNDR)88 should be 17 



instead of 1 as pointed in |^5|, fy in the NDR scheme takes the form 



V 



/ 


3 


-9 






























-9 


3 




































3 


-9 






























-9 


3 




































-1 


3 






























-3 


17 




































-1 


3 






























-3 


17 




































3 


-9 




V 


























-9 


3 


/ 



(42) 



The function Ct-, Cp, and Cg in Eq.(|41|) describe the contributions arising from the 
penguin diagrams of the current- current (5i,2 and the QCD operators Qz-Qq, and the tree- 
level diagram of the magnetic dipole operator Qg, respectively. We here also follow the 
procedure of ref. |]l3l to include the contribution of magnetic gluon penguin operator Qg. 
The effective Wilson coefficients in Eq.(|T) are now renormalization-scheme and -scale 
independent and do not suffer from gauge and infrared problems. The functions Ct-, Cp, 
and Cg are given in the NDR scheme by [|T2], 

"2 

3 

[4 
.3 



C, 



Cr, 



+ ^G(mJ + ^G(m,) 



- - C{mq) - C{mb) 



C3 + 



E 

=u,d,s,c,b 



3 - C{mi 



(C4 + C6), 



(43) 
(44) 



with An 



I + ^G(m„) + ^G(m,)l {C, + 3C2) 

y [O At At 

_ 2m6 

The function G{m, k, fi) is of the form 

2/^2 2 



C{m, k, ji) 



10 2, 
— In — 



9 



+ 



61 

l + 2z, 



(45) 
(46) 

(47) 



^Numerically, such corrections are negligibly small. 

'^The constant term 2/3 in front of C4 + in Cp was missed in ||l^, but recovered firstly in jl^ 
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where z = k"^ /{Aw?), and 




it: 



z <1, 
z> I, 



(48) 



where k"^ is the momentum squared transferred by the gluon, photon or Z to the q'q' pair 
in inclusive three-body decays h qq'q', and m is the mass of internal up-type quark in 
the penguin diagrams. For fc^ > 4m^, an imaginary part of g{z) will appear because of 
the generation of a strong phase at the uu and cc threshold HQ, I 



m 



For the two-body exclusive B meson decays any information on k is lost in the fac- 
torization assumption, and it is not clear what "relevant" fc^ should be taken in nu- 
merical calculation. Based on simple estimates involving two-body kinematics |^9|] or 
the investigations in first paper of ref.[|10|, one usually use the "physical" range for fc^ 



H, H, |T|, |T4| 



~ k 



2 < 



mt 



(49) 



Following refs. |1T2|, |T^, we use k"^ = ml/2 in the numerical calculation and will consider the 
fc^-dependence of branching ratios and CP- violating asymmetries of charmless B meson 
decays. In fact, branching ratios considered here are not sensitive to the value of fc^ within 
the reasonable range of fc^, but the CP- violating asymmetries are sensitive to the variation 
of k\ 



4 Branching ratios oi B ^ PP^ PV decays 

In numerical calculations, we focus on the new physics effects on the branching ratios and 
CP-violating asymmetries for B PP, PV decays. For the standard model part, we will 



follow the procedure of refs.JT2| and compare our SM results with those given in |]T2| , 
Two sets of form factors at the zero momentum transfer from the BSW model |jl5|, as 
well as Lattice QCD and Light-cone QCD sum rules (LQQSR) will be used, respectively. 
Explicit values of these form factors can be found in |[T^ and have also been listed in 
Appendix B. 



Following ||12|, the fifty seven decay channels under study in this paper are also clas- 
sified into five different classes (for more details about classification, see |jl2|) as listed in 
the tables. The first three and last two classes are tree-dominated and penguin-dominated 
decays, respectively. 

• Class-I: including four decay modes, 7r^7r+, p^vr^ and B^ p^K^, the large 
and N^-^^ stable coefficient ai play the major role. 

• Class-II: including ten decay modes, for example B^ tt^tt^, and the relevant 
coefficient for these decays is 02 which shows a strong A^'^-'^-^-dependence. 

• Class-Ill: including nine decay modes involving the interference of class-I and class-II 
decays, such as the decays B^ ttt]'. 
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• Class- IV: including twenty two B PP, PV decay modes such as -B ^ Kr]^'^ 
decays. The amplitudes of these decays involve one (or more) of the dominant 
penguin coefficients 04 6,9 with constructive interference among them. The Class-IV 
decays are N^-^^ stable. 

• Class- V: including twelve B PP, PV decay modes, such as B ^ tt^t]^"' and 
B — »• (pK decays. Since the amplitudes of these decays involve large and delicate 
cancellations due to interference between strong Ai'^-^-^-dependent coefficients 03,5,7,10 
and the dominant penguin coefficients 04,6,9, these decays are generally not stable 
against N^-^^ . 



4.1 Decay amplitudes in BSW model 



With the factorization ansatz |T^, |T2|, , the three-hadron matrix elements or the decay 
amplitude < XY\Heff\B > can be factorized into a sum of products of two current matrix 
elements < X| |0 > and < Y\J2f,\B > ( or < r| Jf |0 > and < X\J2f,\B >). The former 
matrix elements are simply given by the corresponding decay constants fx and .9x1^ 



< 0| J^|X(0-) >= tfxk^, < 0| JJX(I-) >= Mxgxe 



(50) 



where fx {gx ) is the decay constant of pesudoscalar (vector) meson, is the polarization 
vector of the vector meson. For the second matrix element < Y\J2n\B >, the expression 
in terms of Lorentz-scalar form factors |T5|, pO|] , are of the form 

<X(0-)|J^|fi> = 



{ks + kx)^ - 



Ml 



k^ 



^k 



Ml 



A;2 



< xir 



Mb + Mx 
+1 {Mb + Mx 
e* ■ k 



/l/f2 



e*^k1k"xV ~^ {k ) + ie* 



(51) 



f.'^Mx, , ,,,, 



k^ 



A,{k' 



—I 



Mb + M 



X 



{kB + kx) 



Ml -Ml 



k.. 



A^ik' 



(52) 



where k^ = k^ — kx and Mb-, Mx-, My are the masses of meson B, X and Y, respectively. 
The explicit expressions of form factors Fo,i(fc^), Vik"^) and Ao,i,2(/i;^) have been given in 
Appendix B. 

In the generalized factorization ansatz [0, |14|, the effective Wilson coefficients Ci^^ 
will appear in the decay amplitudes in the combinations. 



a2i 



eff 



(Z = l,...,5) 



(53) 



where the effective number of colors N^-I'f is treated as a free parameter varying in the 
range of 2 < N^-l"-^ < 00, in order to get a primary estimation about the size of non- 
factorizable contribution to the hadronic matrix elements. It is evident that the reliability 
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of generalized factorization approach has been improved since the effective Wilson coeffi- 
cients C^^^ appeared in Eq. (|53|) are now gauge invariant and infrared safe. Although N^-l" f 
can in principle vary from channel to channel, but in the energetic two-body hadronic B 
meson decays, it is expected to be process insensitive as supported by the data |]14| . As ar- 
gued in ref. []TB|, N'^^^ {LL) induced by the {V — A){y — A) operators can be rather different 



from Nl^f {LR) generated by {V — A){y + A) operators. Since we here focus on the calcu- 
lation of new physics effects on the studied B meson decays induced by the new penguin 
diagrams in the TC2 model, we will simply assume that N^f^{LL) = N^^^ [LR) = N^^^ 
and consider the variation of N^^^ in the range of 2 < N^^^ < oo. For more details about 
the cases of N^^f{LL) ^ N^^f{LR), one can see for example ref. [0]. We here will also 
not consider the possible effects of final state interaction (FSI) and the contributions from 
annihilation channels although they may play a significant rule for some B PV, VV 
decays. 

The effective coefficients are displayed in the Table |T| and Table ^ for the transitions 
b ^ d {b ^ d ) and b ^ s {b ^ s ), respectively. Theoretical predictions of Oj are made 
by using the input parameters as given in Appendix A and Eq.^, and assuming k"^ = ml/2 
and m^i- = 200GeV. For coefficients 03, ... , aio, the first and second entries in tables 
refer to the values of in the SM and TC2 model respectively. 

The new physics effects on the B decays under study will be included by using the 
modified effective coefficients [i = 3, . . . , 10) as given in the second entries of Table |I| 
and Table ^ In the numerical calculations the input parameters as given in Appendix A, 
B and Eq.@ will be used implicitly. 

From Table |1] and Table 0, one can find several interesting features of coefficients Oj 
because of the inclusion of NP effects: (a) the NP correction to the real part of effective 
coefficients is around 20% for 03,4,5,6, and can be as large as a factor of 4 for coefficients 
07,8,9,10; (b) the NP correction to the imaginary part of is negligibly small; (c) the 
coefficient Oi and 02 remain unchanged since we have neglected the very small tree-level 
NP contributions. 



4.2 Branching ratios of 5 ^ PP decays 

Using above formulas, it is straightforward to find the decay amplitudes of B ^ PP, PV. 
As an example, we present here the decay amplitude M{B' tt^tt^) =< 7i^'n-^\Heff\B^ >, 



M{B- 



with 



G 
"2 

- (a^ + ^(07 - ag) - ^ + (ae - y)^2 ) M, 



(04 + aio + (ae + as)Ri) M^J 



uud 1 



(54) 



Ri 
R2 



2ml 



{mb - mu){mu + m^) 



md{mb - md) ' 



(55) 
(56) 
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Table 1: Numerical values of for the transitions b ^ d \b ^ d ]. The first and second 
entries for 03, ... , Oio refer to the values of in the SM and TC2 model respectively. The 
entries for 03, ... , Oio should be multiplied with 10""^. 







N^ff = 3 


N^ff = 00 






0.995 [0.995] 




1.061 [1.061] 


1.192 [1.192] 


02 




0.201 [0.201] 




0.026 [0.026] 


-0.395 [-0.395] 


as 


-16 


- 7i [-25 - 23z] 




77 [77] 


261 -1- 13« [280 -1- 47^] 




-26 


- 8i [-35 - 24z] 




90 [90] 


322 15i [340 49i] 


04 


-423 - 


-33i [-470 - 117i] 


-467 


-35i [-517- 125i] 


-554 - 39% [-610 - 14h] 




-534- 


- 38i [-580 - 122i] 


-588 


- AOi [-638 - 130i] 


-695 - 45i [-751 - 146i] 


as 


-192 


- 7i [-202 - 23i] 




-71 [-71] 


171 13i [190 47i] 




—195 


— 81 [—205 — 2At\ 




—57 [—57] 


218 15« [237 49«] 




-642 - 


-33i [-689 - 117i] 


-671 


- 35i [-721 - 125i] 


— 798 _ 'iQi r_78zL — IZLI7I 




-718- 


- 38i [-764 - 122i] 


-754 


- 40i [-804 - 130i] 


-827 - 45i [-884 - 146i] 


07 


8.1 - 


-0.9i [7.7- 1.7i] 


6.9 


- 0.9i [6.4 - 1.7i] 


4.3-0.9i [3.9-1.7i] 




34- 


-0.9i [34-1.7i] 


31 


-0.9i [30- 1.7i] 


24.3 - 0.9i [23.9 - 1.7i] 


as 


9.7- 


-0.5i [9.5-0.8i] 


9.0 


- 0.3i [8.8 - 0.6i] 


7.5 [7.5] 




32- 


- 0.5i [31 - 0.8i] 


28 


- 0.3i [27 - 0.6i] 


19.4 [19.4] 


ag 


-83.7- 


-0.9i [-84.1 - 1.7z] 


-90 


- 0.9i [-90 - 1.7i] 


-102-0.9i [-102- 1.7i] 




-153- 


- 0.9i [-153 - 1.7z] 


-164 


- 0.9i [-165 - 1.7i] 


-187-0.9i [-188- 1.7z] 


aio 


-14.4 - 


-0.5i [-14.6 -0.8i] 


-2.6 


-0.3i [-2.5-0.6i] 


37 [37] 




-25 - 


- 0.5i [-25 - 0.8i] 


-6.6 


-0.3i [-6.5-0.6i] 


69 [69] 



Mllf = -t{ml-ml-)UFt{mlo), (57) 



,0 



= -^(m^-m^o)/.Fo«-K-), (58) 

where is the decay constant of vr meson. The form factor F(^'"{m'^) can be found in 
Appendix B. Under the approximation of setting m„ = and m^o = 771^^-, the decay 
amplitude M{B~ tt^tt^) in Eq.(0) will be reduced to the same form as the one given 
in Eq.(80) of [Q: 

M{B- ^ TT-vrO) = -z ^f^FtimlKml - ml) {K.KT. («! + ^2) 

3 1 
-VtbVta X - (-07 + 09 + aio + a8i?2) | (59) 

In the following numerical calculations, we use the decay amplitudes as given in Appendix 
A of ref . [jl2[ directly without further discussions about the details. 

In the B rest frame, the branching ratios of two-body B meson decays can be written 

as 

B{B XY) = ^^\M{B -> XY)\' (60) 
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Table 2: Same as Table |T| but for 6 — > s [6 — * s ] transitions. 





N^ff = 2 






= 3 




= oo 


ai 




0.995 


0.995] 




1.061 


1.061] 


1.192 


[1.192] 


a2 




0.201 


0.201] 




0.026 


0.026] 


-0.395 


[-0.395] 


as 


-21 


- Ui 


-19 - lAi] 




77 


77] 


272 + 29^ 


[269 + 29i\ 




-31 


- 15i 


-30 - 15?] 




90 


90] 


332 + 31i 


[329 + 31i] 




-449 


- 72i 


-442 - 72z] 


-494 


- 77z 


-487 - 77z] 


-585 - 86i 


[-576 - 86z] 




-560 


- 77i 


-553 - 77i] 


-615 


- 82z 


-608 - 82z] 


-725 - 92i 


[-717 - 92z] 


a5 


-198 


- 14i 


-196 - 14i] 




-71 


-71] 


181 + 29i 


[179 + 29i] 




-200 


- 15i 


-199 - 15i] 




-57 


-57] 


229 + 31i 


[226 + 31i] 


ae 


-667 


- 72i 


-661 - 72i] 


-698 


- 77z 


-691 - 77i] 


-758 - 86i 


[-750 - 86i] 




—744 


— 77i 


—737 — 77^] 


-782 


— 82i 


-774 - %2i\ 


-858 - 92i 


[-850 - 92^] 


dj 


7.9 - 


- 1.3i 


7.9- IM] 


P. P. 

0.0 


— i.Oi 


0. ( — i.OiJ 


/II 1 

41. i — l.Ot 


'A 9 1 

41. Z — i.OiJ 




34- 


- 1.3i 


34 - 1.3i] 


31 


- \M 


31 - 1.3i] 


24- 1.3i 


[24 - 1.3i] 


as 


9.6 - 


- 0.6i 


9.6 - 0.6i] 


8.9 


-0.4i 


8.9 - 0.4i] 


7.5 


[7.5] 




32- 


- 0.6i 


32 - 0.6i] 


28 


- 0.4i 


28 - 0.42] 


19.4 


[19.4] 


ag 


-84- 


- 1.3z 


-84-1.3i] 


-90 


- \M 


-90- 1.3i] 


-102 - 1.32 


-102 - 1.3i] 




-153- 


- IM 


-153- 1.3i] 


-165 


- \M 


-164- 1.3i] 


-188 - 1.3i 


-187- 1.3i] 




-14.5 - 


- 0.6i 


-14.5 - 0.6i] 


-2.2 


-0.4z 


-2.6 - 0.4i] 


37 


[37] 




-25- 


- 0.6i 


-25 - 0.6i] 


-6.6 


-0.4z 


-6.6 - 0.4z] 


69 


[69] 



for B PP decays, and 

B{B XY) = ^^mB XY)/{e ■ ps)\' (61) 
1 tot 8ttM^ 

for B ^ PV decays. Here TtotiB-) = 3.989 x lO^^^ Q^y ^nd 1^4(5°) = 4.219 x lO-^^^GeV 
obtained by using t{B~) = 1.65j9s and t{B'^) = 1.56ps pb is the four-momentum of 



the B meson, My and e is the mass and polarization vector of the produced light vector 
meson respectively, and 

1 



bl = ^ylM^B - (Mx + My)2][M| - (Mx - My)2] (62) 

is the magnitude of momentum of particle X and Y in the B rest frame. 

In Tables 0-^, we present the numerical results of the branching ratios for the twenty 
B PP decays and thirty seven B PV decays in the framework of the SM and 
TC2 model. The theoretical predictions are made by using the central values of input 
parameters as given in Eq.(|]) and Appendix A and B, and assuming m^f = 200GeV 
and N^-^'^ = 2, 3, oo in the generalized factorization approach. The fc^-dependence of the 
branching ratios is weak in the range of fc^ = ml/2 ± 2 GeV^ and hence the numerical 
results are given by fixing k"^ = ml/2. The currently available CLEO data|]19|, ^ ^ are 
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listed in the last column. The branching ratios collected in the tables are the averages 
of the branching ratios of B and anti-i? decays. The ratio 6B describes the new physics 
corrections on the SM predictions of corresponding branching ratios and is defined as 

6B{B XY) = ^ XY)SM • (63) 

By comparing the numerical results with the CLEO data, the following general features 
of i? ^ PP decays can be understood: 

• The SM predictions for five measured B^ tt+tt^ and B Kn decay modes 
are consistent with the CLEO data. But for the measured B Krj' decays, the 
observed branching ratio are clearly much larger than the SM predictions |TT], |T2], Q . 
All other estimated branching ratios in Table ^ and Table ^ are consistent with the 
new CLEO upper limits. 

• The uncertainties of the SM predictions for the branching ratios oi B ^ PP decays 
induced by varying k"^ is ~ 10% within the range of k"^ = ml/2 ± 2GeV^. 

• For most class-II, IV and V decay channels, such as B ^ r]r]^'\ B Ktt,B Kr]', 
etc, the NP enhancements to the decay rates can be rather large: from 20% to 70% 
w.r.t the SM predictions. 

• For most B PP decay channels, the magnitude of NP effects is insensitive to the 
variations of m^;- and N^^^ . 

• The central values of the branching ratios obtained by using the LQQSR form factors 
will be generally increased by about 15% when compared with the results using the 
BSW form factors, as can be seen from Table | and Table |. No matter the BSW 
or the LQQSR form factors was used, the magnitude and whole pattern of the new 
physics corrections to the decay rates in study remain basically unchanged. 

• Both new gluonic and electroweak penguin diagrams contribute effectively to most 
decay modes. 



4.2.1 B TTTY, Kn decays 

There are so far seven measured B —>■ PP decay modes EO, |21], p5|: 



B(B^n-^rr-) - / (4-31^:1 ± 0.5) x lO"*^ [CLEO], 

BiB^K^n'^) - / (ll-6;f7;i:^)xl0-6 [CLEO], 
B{B^K n) - j (18.8+5.5 ±2.3) X 10-6 [Belle], 

r (17.2i|t± 1.2) X 10-6 [CLEO], 
B{B^K+n-) = I (12.5i|S] il:?± 2.3) X 10-6 [BaBar], (66) 

[ (17.4t|^ ± 3.4) X 10-6 [Belle], 

B{B^K%+) = (18.2+^6 ± 1.6) X 10-6 [CLEO], (67) 
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Table 3: i? — > PP branching ratios (in units of 10 ®) using the BSW form factors, with 
P = ml/2, A = 0.81, A = 0.2205, p = 0.12, rj = 0.34, N^^f = 2, 3, oo and assuming = 
200 GeV, in the SM and TC2 model by employing generalized factorization approach. The 
last column contains measured branching ratios and upper limits (90%C.L.) [19,20] 









SM 


TC2 


5B [%] 




Channel 


Ulass 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


Data 






T 

I 


9.10 


10.3 


13.0 


9.27 


10.5 


13.2 


1.9 


1.8 


1.6 


4.3_i;5 =t 0.5 


■r-»n 




TT 

11 


0.28 


0.15 


0.92 


0.28 


0.16 


0.94 


1.0 


6.3 


2.8 


< 9.3 


B+ - 


- 7r+7r° 


TTT 

111 


6.41 


5.06 


2.85 


6.41 


5.07 


2.85 


0.1 


0.1 


0.1 


< 12.7 




rjf] 


TT 
ii 


0.14 


0.10 


0.29 


0.20 


0.17 


0.38 


40 


64 


30 


< lo 


5° 


rji]' 


TT 

11 


0.14 


0.08 


0.38 


0.19 


0.13 


0.45 


30 


67 


19 


< 27 




rj'rj' 


TT 

11 


0.05 


0.01 


0.13 


0.04 


0.02 


0.14 


13 


73 


7.8 


< 47 


B+ - 


■ TT+ry 


III 


3.51 


2.78 


1.75 


3.85 


3.17 


2.25 


10 


14 


28 


< 5.7 


Jd — 


^+^' 

- 711] 


TTT 
ill 


2.49 


1.88 


1.02 


2.59 


1.99 


1.16 


3.8 


5.8 


13 


^ 1 o 

< IZ 


B^^ 




V 


U.ZD 


n on 

u.zy 


n on 


U.OO 


U.4z 


U.O i 




A A 

44 


A C 

4o 


< 2.9 


B^^ 


n\ 


V 


0.06 


0.08 


0.14 


0.08 


0.10 


0.18 


37 


35 


26 


< 5.7 


B+ - 




IV 


12.0 


13.5 


16.7 


19.6 


21.8 


26.5 


63 


61 


59 


-L J- -0-2.7-1.3 


B^^ 




IV 


17.8 


19.8 


24.0 


24.4 


26.9 


32.2 


37 


36 


35 


17.2l|t± 1.2 


B+ - 




IV 


19.9 


23.2 


30.6 


27.7 


32.7 


44.0 


39 


41 


44 


18.2+to ± 1-6 






IV 


7.27 


8.31 


10.7 


7.95 


9.36 


12.6 


9.3 


13 


18 




B+ - 


► K+T] 


IV 


3.91 


4.56 


6.07 


4.09 


5.08 


7.45 


4.6 


11 


23 


< 6.9 


B+ - 


' K+ri' 


IV 


22.6 


28.5 


42.4 


33.8 


41.6 


59.5 


50 


46 


40 


80l^° ± 7 




K^-q 


IV 


3.22 


3.63 


4.58 


3.33 


3.90 


5.23 


3.6 


7.5 


14 


< 9.3 




K°ri' 


IV 


21.9 


28.2 


43.0 


32.9 


41.3 


61.2 


50 


47 


43 


89lli ± 9 


B+ - 




IV 


1.16 


1.35 


1.78 


1.61 


1.90 


2.55 


38 


40 


43 


< 5.1 


B^^ 




IV 


1.10 


1.28 


1.68 


1.52 


1.80 


2.41 


38 


40 


43 


< 17 



n(B^K'^'\ - / (14-6i|?tli)xl0-6 [CLEO], 

B{B^Kn) - j (21+9-3 ^ 10-6 p^He], ^^^^ 

l^[B^K T)) - I (62 ±18 ±8) X 10-6 [BaBar], 

BiB^K'^T]') = (89l{^ ± 9) X 10-6 [CLEO], (70) 

The measurements of CLEO, BaBar and Belle are in good agreement within errors. 

As a Class-I decay channel, the B^ — >■ 7r"'"7r- decay are dominated by the b u tree 
diagram. This mode together with B^ tt'^tt^ and B~^ tt+tt^ decays play an important 
role in determination of angle a. For all three 5 — > tttt decay modes, the new penguin 
enhancement is very small, < 6.3% for N^^f = 2 — oo, as hsted in tables 3 and 4. The 
theoretical predictions in the SM and TC2 model are consistent with the CLEO data. 

For B^ Tj^'^T)^'^ decays, the NP enhancement is varying in the range of 10% to 
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Table 4: Same as Table ^ but using the LQQSR form factors. 









SM 


TC2 


5B [%] 




Channel 


olass 


2 


3 


00 


2 


3 


00 


2 


3 


00 


Jjata 


- 




T 
i 


10.8 


12.3 


15.5 


11.0 


12.5 


15.8 


1.9 


1.8 


1.6 


A + 1. b 1 n rr 


- 




TT 
ii 


0.33 


0.18 


1.09 


0.33 


0.19 


1.12 


1.0 


6.3 


2.8 


^ n 


B+ - 


7r+7r° 


TTT 
ill 


7.62 


6.02 


3.39 


7.63 


6.03 


3.39 


0.1 


0.1 


0.1 


^19 7 


5° - 


7777 


TT 
ii 


0.17 


0.13 


0.36 


0.24 


0.21 


0.47 


40 


64 


30 


< lo 


B^ - 


rjri' 


TT 
ii 


0.17 


0.09 


0.45 


0.22 


0.15 


0.53 


30 


67 


19 


< l( 


5° - 




TT 
ii 


0.05 


0.01 


0.15 


0.05 


0.02 


0.16 


13 


73 


7.8 


< 47 


B+ - 




TTT 
ill 


4.25 


3.37 


2.13 


4.66 


3.83 


2.73 


9.6 


14 


28 


< 0. ( 


B+ - 


TT + T]' 


TTT 
iii 


2.90 


2.17 


1.17 


3.01 


2.30 


1.33 


3.8 


5.8 


14 


^19 


B" - 


n 


V 


0.31 


0.35 


0.47 


0.43 


0.50 


0.69 


42 


44 


46 


^ 9 n 

< z.y 


S° - 




V 


0.07 


0.09 


0.17 


0.09 


0.12 


0.21 


37 


36 


26 


< 0. / 


£3 


-> A TT 


T\/ 
1 V 


14.3 


16.0 


19.8 


23.2 


25.8 


31.4 


63 


61 


58 


11 a+'SM+lA 
ii.O_2.7_1. 3 






TA T 
iV 


Zi.Z 


ZO.O 


Zo.O 


9Q n 


Q9 n 


00.4 


Q7 
( 


00 


oO 


1 ^7 + 2.5 _1_ 1 

1 ( .2_2.4 ± i.Z 


B+ - 




IV 


23.7 


27.7 


36.4 


33.0 


38.9 


52.3 


39 


41 


43 


18.2+40 ± 1.6 


fiO - 




IV 


8.68 


9.92 


12.7 


9.51 


11.2 


15.1 


9.6 


13 


18 


14.611? tlj 


fi+ - 


K+r] 


IV 


4.37 


5.10 


6.80 


4.54 


5.66 


8.33 


3.9 


11 


22 


< 6.9 


B+ - 


K+T]' 


IV 


26.2 


33.1 


49.2 


39.2 


48.2 


69.1 


50 


46 


40 


80tl^ ± 7 


B' - 


K^T] 


IV 


3.57 


4.02 


5.07 


3.67 


4.30 


5.76 


2.8 


6.8 


14 


< 9.3 


5°- 


^ K'>7]' 


IV 


25.5 


32.7 


49.9 


38.1 


48.0 


71.1 


50 


47 


43 


891}^ ± 9 


S+ - 




IV 


1.35 


1.58 


2.07 


1.87 


2.21 


2.96 


38 


40 


43 


< 5.1 






IV 


1.28 


1.49 


1.96 


1.77 


2.09 


2.80 


38 


40 


43 


< 17 



70%. For B^ -n^T]^'^ decays, the NP enhancement is around 10% and depend on N^^^ 
moderately. For 5° tt^t]^ ^ decays, the NP enhancement is large, 30% — 60%, and 
insensitive to the variation of N^^f . 

In the SM, the four Class-IV decays B — > Kt^ are dominated by the b ^ sg gluonic 
penguin diagram, with additional contributions from b ^ u tree and electroweak penguin 
diagrams. Measurements of i? — > Ktc decays are particularly important to measure the 
angle 7. In the TC2 model, the new penguin diagrams will interfere with their SM 
counterparts and change the SM predictions for the branching ratios and CP-violating 
asymmetries. 

It is well known that the effective Hamiltonian calculations of charmless hadronic 
B meson decays contain many uncertainties including form factors, light quark masses, 
CKM matrix elements, QCD scale and external momentum fc^. As a simple illustration of 
the theoretical uncertainties, we calculate and show the branching ratios of four B — * Kir 
decay modes by using F^'^{0) = 0.20 ( preferred by the CLEO measurement of i? ^ tt+tt" 
mode fSll ) instead of the ordinary BSW value F^'"{0) = 0.33 ( all other input parameters 
remain unchanged ) and by varying rj, k"^ and m^p in the ranges oi r] = 0.34 ± 0.08, 
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Table 5: Same as Table H, but for branching ratios of 5 ^ PP decays with new physics 
contributions from charged-scalar gluonic penguins only. 









TC2: 


QCD only 


6B [%] 






Channel 


Class 


2 


3 


oo 


2 


3 


oo 


Data 


5° 


— > TT + TT 


1 


9.27 


10.5 


13.3 


1.90 


1.86 


1.90 


4.3li:^±o.5 


5° 




11 


0.34 


0.22 


1.01 


23.3 


47.9 


9.85 


< 9.3 


B+ 




TTT 

111 


6.41 


5.06 


2.85 


0.0 


0.0 


0.0 


< 12.7 




rji] 


11 


0.17 


0.14 


0.34 


23.8 


37.1 


16.5 


< 18 




rjrj' 


11 


0.19 


0.13 


0.45 


32.1 


68.7 


17.7 


< 27 


fio 


7]' 7]' 


11 


0.05 


0.02 


0.15 


29.3 


162 


16.9 


< 47 


B+ 


— s> Ti^rj 


111 


3.75 


3.05 


2.10 


6.86 


9.86 


19.7 


< 5.7 


B+ 


— s> TT+r]' 


111 


2.65 


2.05 


1.25 


6.11 


9.38 


21.8 


< 12 


B° 


— s> Ti^rj 


V 


0.36 


0.41 


0.54 


40.2 


40.5 


37.8 


< 2.9 






V 


0.12 


0.15 


0.24 


107 


97.2 


65.2 


< 5.7 


B^ 


— > K^-K^' 


IV 


16.0 


18.0 


22.4 


33.0 


33.5 


34.1 


1 1 n + 'i 0+1 4 
11-612:7^1:3 


5° 




IV 


24.5 


27.3 


33.3 


37.7 


38.2 


39.1 


17.2t|| ± 1.2 


B+ 




1 V 


97 


31.7 


41.5 


37.0 


36.5 


oo. ( 


io.z 4 Q in: i.u 


5° 




IV 


10.4 


11.8 


15.1 


42.5 


42.4 


42.0 


i4.6ii:?;it 


5+ 


K+T] 


IV 


5.02 


5.90 


7.95 


28.4 


29.4 


30.8 


< 6.9 


B+ 


K+7]' 


IV 


37.4 


45.2 


63.2 


65.6 


58.8 


49.0 


80tl^ ± 7 


B' 


K^T] 


IV 


4.25 


4.85 


6.22 


32.1 


33.6 


35.8 


< 9.3 


5° 




IV 


36.1 


44.2 


63.1 


64.4 


57.1 


46.9 


89t\l ± 9 


B+ 




IV 


1.59 


1.84 


2.41 


36.4 


35.9 


35.2 


< 5.1 


5° 




IV 


1.50 


1.74 


2.27 


36.4 


35.9 


35.2 


< 17 



A;2 = ml/2 ± 2GeV'^, = 200 ± 100 GeV, and setting N^^f = 2, 3, 00: 

/ (5.8 ±0.1^^it^;^7)x 10-6 in SM, 

" \ (10.1 ±o.iiJ:° X 10-6 in tc2, ^^^^ 

j (7.3±0.li°:|lo:i)xlO-6 in SM, 
= (9.9 ±0.llJ:^lJ:rjL)x 10-6 in TC2, ^^'^ 



B{B - 




B{B - 




B{B - 




B{B - 





(8.5±0.0i?|lf^) X 10-6 in SM, 

(12.0 ±0.0^1:^ its -as) X 10"' in TC2, 
(2.5 ± O.Olgi Ig-^) X 10-6 in SM, 
(2.3±0.0l°:^l°|;°i)xl0-6 in TC2, 



(73) 
(74) 



where the first, second and third error correspond to the uncertainty 6ri = ±0.08, 6q^ = ±2 
and 2 < N^-^-^ < 00 respectively, while the fourth error refers to 5m^i- = ±100 GeV. By 
comparing the ratios in tables (|^, ^) and in Eqs. ([7ll - [74D , it is easy to see that the central 
values of the branching ratios B{B — > Kn) are greatly reduced by using Fq^'^(O) = 0.20 
instead of 0.33, the new physics enhancements therefore become essential to make the 
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theoretical predictions being consistent with data. 

Fig.|] shows the mass and A^^-^-^-dependence of the ratios B{B K^tt'^ ) in the SM 
and TC2 model using the input parameters as given in Appendix A and B and employing 
the BSW form factors. In Fig. (2a), we set N^'^-^ = 2 and assume that = 100 — SOOGeV. 
In Fig. (2b), we set = 200GeV and assume that ^ = l/N^^f = - 0.5. In Fig.(g) the 
short-dashed line and solid curve show the branching ratio of K^ix^ decay in the 

SM and TC2 model, respectively. The band between two dots lines corresponds to the 
CLEO data with 2a errors: B{B K+n^) = {U.Qtt^o) x 10^^. 

Same as Fig.^ Figs. (3,4,5) show the mass and A"^-^ ■^-dependence of the branching 
ratios of decay B — >■ K~^7r^ , K^7r~ , K^7r~^ and K^n^, respectively. In these three figures, 
the short-dashed lines and solid curves show the branching ratios for relevant decay modes 
in the SM and TC2 model. The band again refers to the corresponding CLEO data with 
2a errors, respectively: B{B K+n-) = (17.2l|^) x 10-^ B{B K%+) = (18.2t|^) x 
lO"'' and B{B K^ir^) = {I4:.6tlli) x 10"*^. The large theoretical uncertainties are not 
shown in all four figures. 

Although the new physics enhancements to the branching ratios of -B — K^tt and 
K^7r~^ decays are relatively large as illustrated in Figs. (2,3,4), the theoretical predictions 
for B Kit decays in the TC2 model are still consistent with the CLEO measurements 
within 2a errors after taking into account existed large theoretical uncertainties. If one 
uses F^'^{0) ~ 0.20 instead of 0.33, the new physics effects will play an important role to 
boost the theoretical predictions for branching ratios of i? ^ Kir decays. 



4.2.2 B Kt]'''^ decays and new physics effects 

In the SM, the Class-IV decays B Kr]^'^ are expected to proceed primarily through 
b ^ s penguin diagrams and h u tree diagram. In TC2 model, the new gluonic and 
electroweak penguins will also contribute through interference with their SM counterparts. 
The CLEO data of 5 — >■ Krf^ ^ decays with recent measurements of 5 — >■ tttt, fCvr, etc, 
provide important constraints on the theoretical picture for these charmless B meson 
decays. 

For B^ K^rj and B'^ K^rj decay modes, the new physics enhancement is less 
than 10% for N^^^ ~ 3. The theoretical predictions in both the SM and TC2 model are 
consistent with the new CLEO upper limits: B{B — > K^rf) < 6.9 x 10~^ and B{B 
K^r]) < 9.3 X lO-^ig. 

For B Kr]' decay modes, the situation is very interesting now. Unexpectedly 



large B Kr]' rates were firstly reported by CLEO in 1997|5^, and confirmed very 
recently pO|, The Krj' signal is large, stable and has small errors (~ 14%). Those 
measured ratios as given in Eqs.(^|7D|) are clearly much larger than the SM predictions 
(the contribution from the decay h s(cc) s{ri,ri') have been included [0, [1^ ) as 
given in tables (UM and illustrated by the short-dashed line in Figs.(P,IW) where only the 



central values of theoretical predictions are shown. Furthermore, Lipkin's sum rule 



B{K+ri') + B{K+ri) = B{K+7r^) + i3(i^°7r+) (75) 



is also strongly violated (~ Aa) [|^: 82.2^11'^ = 29.8lll'l. At present, it is indeed difficult 



to explain the observed large rate for B Krj' in the framework of SM [^O], This 
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fact strongly suggest the requirement for additional contributions unique to the r/' meson 

ml/2 ± 2Gey2, 



in the framework of the SM, or from new physics beyond the SM |^ 
By varying r/, fc^ and in the ranges of r/ = 0.34 ± 0.08, /c^ 
= 200 ± 100 GeV, and setting N^f^ 



2, 3, oo, we find that 



B{B K+7]') 



+2.7 +13.9-^ 
-2.2 -6.9 



6.2 +17.9 +3.3\ 

4.3 -7.8 -2.7 ) 



(26.5 ±0.1 
(41.6 ±0.1 

(28.2 ±0.1^111^^3 
(41.3 ±0.1111+^^-9+=^ 



X 10- 



X 10 



X 10" 

6 



2.7 



X lo- 



in SM, 

in TC2, 

in SM, 

in TC2, 



(76) 
(77) 



where the first to the fourth error corresponds to the uncertainty drj = ±0.08, = ±2 
and 2 < N^^^ < oo and 5m^f = ±100 GeV, respectively. If we use the LQQSR form factors 
instead of the BSW form factors, the central values of BR{B — > Kr]^'^) will be increased 
by about 15%. The NP enhancements to B ^ Krj' decays are significant numerically, 
~ 50% for = 200 GeV. 

Taking into account all uncertainties considered here, the theoretical predictions for 
the magnitude of B(i? Krj') in the SM and TC2 model are 



B{B K+7]') 
B{B K%') 



(20- 


53) 


X 


10^^ 


in 


SM, 


(30- 


74) 


X 


10-6 


in 


TC2 


(19- 


52) 


X 


10-^ 


in 


SM, 


(28- 


76) 


X 




in 


TC2 



(78) 
(79) 



It is evident that the theoretical predictions for ratios B{B Kr]') become now consistent 
with the CLEO data due to the NP enhancements. This is a plausible new physics 
interpretation for the large B — »• Krj' decay rates. 

Figs.(P,|^) show the mass and N^^^'-f^ dependence of the ratios B{B — >■ Kt]') in the SM 
and TC2 model using the input parameters as given in Appendix A and B and employing 
the BSW form factors. The short-dashed and solid curves in Figs. (6,7) show the central 
values of theoretical predictions. The band corresponds to the CLEO measurements with 
2a errors. 



4.3 Branching ratios of B ^ PV decays 



In tables (|g-|g) we present the branching ratios for the thirty seven B PV decay modes 
involving b ^ d and b s transitions in the SM and TC2 model by using the BSW and 
LQQSR form factors and by employing generalized factorization approach. Theoretical 
predictions are made by using the same input parameters as those for the B PP decays 
in last subsection. The measured branching ratios from CLEO [|19|, |2^ for six B PV 
decay modes, B p+vr^, p°7r+, ujtt^, K*^r], K*^ri, K*^^^ , have been given in last column 
of Table |^. BaBar and Belle also reported their measurements for B{B^ 



p vr^ 



and B{B 



B{B^ - 
B{B^ 



(49 ± X 10-6 

(27.6l?:H4.2) 



[BaBar] , 
X 10-6 [CLEO], 

(17.2l|^± 1.8) X 10-6 [Belle]. 



(80) 
(81) 
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Table 6: -B — > PV branching ratios (in units of 10 ^) using the BSW [ LQQSR ] form 
factors in the SM. The last column shows the CLEO and Belle measurements and upper 
limits (90% C.L.) [19,20,25]. 



N^ff = 2 
21.1 [25.1J 
6.5] 
0.58] 
6.63] 
16.0] 



N^^f = oo 
30.3 [36.0] 
8.19 [9.4] 
2.05 [2.41] 
0.71 [0.78] 
10.9 [13.1] 



Channel 



Class 



= 3 
283f 
[7.4] 
0.07] 
3.97] 
15.0] 



Data 



5° - 
B+ 



P ' TT 



I 
I 

II 

III 

III 



5.7 
0.49 

5.72 
13.5 



24.0 
6.48 
0.06 
3.46 
12.6 



} 27.6I?:4 ± 4.2 
< 5.1 

10.4l|^±2.1 

<43 



- 

B°- 

B+ 

B+ 



p^'i 

► p^T] 

> p+ri' 



II 
II 
III 
III 



0.01 
0.02 
5.44 
4.35 



0.02] 
0.01] 
6.57] 
5.02] 



0.02 
0.002 
4.75 
3.81 



0.02] 
0.003] 
5.79] 
4.40] 



0.06 
0.03 
3.54 
2.85 



0.08] 
0.03] 
4.38] 
3.29] 



< 10 

< 12 

< 15 

< 33 



B^ ^UTT 
B+ UJTl 

B° 
5° 







our) 



II 
III 
II 
II 



0.29 
6.32 
0.32 
0.20 



0.35] 
7.35] 
0.38] 
0.23] 



0.08 
3.75 
0.03 
0.001 



0.09] 
4.31] 
0.04] 
0.002] 



0.15 
0.78 
0.82 
0.68 



0.19] 
0.85] 
0.98] 
0.79] 



11.3] 



< 5.5 

'+2:9 ±1-4 

< 12 

< 60 



50 ^ 07rO 
B+ 07r+ 

B^ 4)7]' 



V 
V 
V 
V 



0.03 
0.06 
0.01 
0.01 



0.04] 
0.08] 
0.01] 
0.01] 



0.002 
0.004 
0.001 
0.001 



0.002] 
0.005] 
0.001] 
0.001] 



0.23 
0.49 
0.09 
0.07 



0.27] 
0.58] 
0.10] 
0.081 



< 5.4 

< 4 

< 9 

< 31 



B+ K*^K+ 

B+ K*+K^> 
B° ^ K*^K^ 



IV 

IV 

V 

IV 



0.42 
0.40 
0.004 
0.004 



0.49] 
0.46] 
0.006] 
0.006] 



0.53 
0.50 
0.002 
0.002 



0.61] 
0.58] 
0.003] 
0.003] 



0.78 
0.73 
0.001 
0.001 



0.90] 
0.85] 
0.001] 
0.001] 



< 5.3 



< 12 



50 ^ p^K^ 
B+ ^ p°K+ 
B° p-K+ 
B+ ^ p+i^o 



IV 
IV 
I 

IV 



0.52 
0.39 
0.54 
0.11 



0.60] 
0.46] 
0.62] 
0.12] 



0.53 
0.31 
0.59 
0.05 



0.62] 
0.36] 
0.68] 
0.05] 



0.72 
0.31 
0.70 
0.005 



0.83] 
0.36] 
0.81] 
0.006] 



< 27 

< 17 

< 25 

< 48 



26.41^;^ ± 3.3 

< 35 

13.8iy ± 1.6 

< 20 



B+ K*+r] 
B+ K*+r]' 
B° K*^rj 



IV 
III 
IV 
V 



2.43 
0.66 
2.32 
0.33 



3.12] 
1.14] 
2.98] 
0.69] 



2.39 
0.36 
2.54 
0.09 



3.04] 
0.61] 
3.23] 
0.23] 



2.32 
0.24 
3.06 
0.31 



2.89] 
0.23] 
3.82] 
0.26] 



+8+4 
6-5 



50 ^ K*+7r- 

B+ K*+7r° 
B+ K*^n+ 



IV 
IV 
IV 
IV 



8.59 
2.44 
4.95 
7.35 



10.2] 
2.77] 
6.09] 
8.75] 



9.67 
3.02 
5.55 
9.23 



11.5] 
3.43] 
6.84] 
11.0] 



12.0 
4.42 
6.91 
13.6 



14.3] 
5.01] 
8.52] 
16.2] 
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< 3.6 

< 31 

< 16 



< 17.21^:^ ± 1.8 



B+ 

B^- 



> (l)K+ 



V 

V 



22.1 

20.9 



25.7] 

24.3] 



11.5 

10.9 



13.4] 

12.6] 



0.60 

0.57 



0.70] 

0.66] 



< 2i 



B+ 



V 
V 



3.31 
3.53 



3.86] 
4.11] 



0.002 
0.25 



0.003] 
0.28] 



13.3 
16.5 



15.4] 
19.2] 



< 21 

< 7.9 
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Table 7: B ^ PV branching ratios 
factors in the TC2 model. 



(in units of IQ-^) using the BSW [LQQSR] form 



TC2 



SB \% 



Channel 



Class 



B+ 



oo 



oo 



+ 







1 
I 

II 

III 
111 



21.2 
5.70 
0.49 

5.72 
13.6 



25.3] 
'6.54] 
0.58] 
6.63] 
16.2] 



24.1 
6.48 
0.06 

3.46 
12.7 



28.7] 
7.44] 
0.07] 

3.98] 
15.1] 



30.4 
8.19 
2.06 

0.73 
11.0 



36.2 
9.40 
2.43^ 
0.81 
13.2 



0.71 
-0.06 
0.05 

-0.02 
0.83 



0.63 
-0.05 
5.90 
0.19 
0.86 



0.50 
-0.03 
0.53 
3.62 
0.92 



B°- 
B+ 

B+ 



B+ 
B° 
5° 



p'^i 



11 
11 
III 
111 



0.03 
0.01 
5.49 
4.35 



IT 



0.04] 
0.01] 
6.63] 
5.01] 



0.03 
0.002 
4.82 
3.81 



0.04] 
0.002] 
5.86] 
4.39] 



0.09 
0.03 
3.62 
2.85 



0.11 
0.03 
4.48 
3.29 



96.7 
-21.5 

0.96 
-0.19 



116 
-30.5 

1.29 
-0.08 



49.1 
2.39 
2.26 
0.15 



LOTY 
> LUTT 

CUT)' 



11 
111 
11 
11 



0.33 
6.58 
0.38 
0.21 



0.39] 
7.65] 
0.45] 
0.24] 



0.08 
3.84 
0.06 
0.002 



0.10] 
4.43] 
0.07] 
0.002] 



0.17 
0.78 
0.82 
0.69 



0.21 
0.85 
0.98 
0.79 



12.3 
4.05 
19.0 
6.00 



3.23 
2.63 
107 
63.8 



13.3 
-0.10 
-0.43 

0.94 



B°- 

B+ 

5°- 



V 
V 
V 
V 



0.03 
0.06 
0.01 
0.01 



B° K*°K° 

B+ K*+K^ 
5° K*°K° 



0.03] 
0.07] 
0.01] 
0.01] 



0.009 
0.02 
0.003 
0.003 



[0.01] 
0.02] 
0.004] 
0.003] 



0.37 
0.79 
0.14 
0.11 



0.44 
0.94^ 
0.17 
0.13 



-9.1 
-9.1 
-9.1 
-9.1 



351 

351 
351 
351 



62.8 

62.8 
62.8 
62.8 



IV 
IV 

V 
IV 



0.64 
0.60 

0.00 
0.001 



0.74] 
0.70] 

[0.002] 
0.002] 



0.81 
0.77 

0.005 
0.005 



0.95] 
0.89] 

0.001] 
0.001] 



1.22 
1.16 

0.01 
0.01 



1.43 
1.35 

0.01 
0.01 



52.1 
52.1 

-71.S 
-71.S 



54.4 
54.4 

-72.9 
-72.9 



57.7 
57.7 

886 
885 



5° ^ p°K° 
B+ ^ 

5° ^ p-K+ 
B+ ^ p+i^o 



IV 
IV 
I 

IV 



0.85 
0.86 
0.38 
0.04 



0.99] 
1.00] 
0.44] 
0.04] 



0.92 
0.92 
0.45 
0.005 



1.07] 
1.07] 
0.51] 
0.006] 



1.21 
1.24 
0.60 

0.09 



1.41 
1.44; 
0.69 
0.11 



64.5 
118 
-29.6 
-66.0 



72.0 
198 
-24.5 
-89.7 



69.3 
299 
-14.3 

1686. 



B+ K*+7] 
B+ K*+ri' 

50 ^ ^*o^/ 



IV 
111 
IV 
V 



4.02 
0.32 
3.74 
0.10 



5.18] 
0.55] 
4.82] 
0.24] 



3.81 
0.24 
4.07 
0.10 



4.85] 
0.31] 
5.18] 
0.08] 



3.39 
0.50 
4.79 
0.93 



4.23 
0.44 
5.97 
0.87^ 



65.4 
-50.8 

61.4 
-71.3 



59.0 
-32.8 
59.8 
2.18 



46.1 
108 
56.6 
196 



B" 
B^ 
B+ 
B+ 



K*+n- 



IV 
IV 
IV 
IV 



13.6 
2.74 
9.38 
11.2 



16.2] 
2.97] 
11.7] 
13.4] 



14.7 
3.58 
10.2 
14.3 



17.5] 
3.89] 
12.8] 
17.1] 



17.1 
5.63 
12.0 
21.6 



20.4 
6.17 
15.1 
25.7 



58.1 
12.5 
89.4 
53.0 



52.3 
18.5 
84.0 
55.2 



42.7 
27.4 
74.3 
58.5 



B+ 
B^- 



► (f)K+ 



V 

V 



29.4 

27.8 



34.3] 
32.4] 



15.3 
14.5 



17.9] 
16.9] 



0.82 
0.77 



0.95 
0.90 



33.3 
33.3 



33.5 
33.5 



35.2 
35.2 



B+ 



>uK+ 



V 
V 



4.49 
5.18 



5.23] 
6.04] 



0.003 
0.24 



0.003] 
0.27] 



18.0 
22.8 



21.0 
26.5 



35.6 
46.9 



12.7 
-4.63 



35.9 
38.3 
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Table 8: S — > PV branching ratios (in units of 10 ^) using the BSW form factors in TC2 
model with new contributions induced by charged-Higgs gluonic penguin diagrams only. 









TC2 


: QCD only 


5B [%] 


Channel 


Class 


2 


3 


oo 


2 


3 


oo 




r 


I 


21.2 


24.1 


30.4 


0.64 


0.62 


0.59 


B^^ 


r 


I 


5.70 


6.48 


8.19 


-0.06 


-0.05 


-0.04 


5°^ 




II 


0.54 


0.11 


2.12 


9.68 


95.0 


3.12 


B+ - 


r 


III 


5.77 


3.52 


0.78 


0.94 


1.75 


10.6 


B+ - 


r 


III 


13.6 


12.7 


11.0 


0.32 


0.38 


0.54 


B» ^ 


r 1 


II 


0.03 


0.03 


0.08 


105 


115 


41.1 


B^^ 




II 


0.004 


0.003 


0.04 


-51.1 


-9.35 


25.5 


B+ - 


1 1 


III 


5.47 


4.79 


3.59 


0.61 


0.82 


1.44 


B+ - 


r 1 


III 


4.34 


3.81 


2.86 


-0.21 


0.0 


0.56 


5° ^ 


UJTT^ 


II 


0.43 


0.13 


0.15 


46.8 


70.6 


-1.97 


B+ - 




III 


6.59 


3.85 


0.77 


4.33 


2.87 


-0.41 


B^^ 


our] 


II 


0.37 


0.05 


0.82 


16.0 


73.8 


-0.08 


B^^ 


ujr]' 


II 


0.23 


0.006 


0.68 


13.6 


363 


-0.89 




T 


V 


0.04 


0.002 


0.30 


39.6 


13.5 


32.2 






V 


0.09 


0.005 


0.64 


39.6 


13.5 


32.2 


B°^ 


(brj 


V 


0.02 


0.001 


0.11 


39.6 


13.5 


32.2 


I?" ^ 




V 


0.01 


0.001 


0.09 


39.6 


13.5 


32.2 


B+ - 




IV 


0.65 


0.79 


1.13 


54.9 


51.0 


45.2 


B^^ 




IV 


0.61 


0.75 


1.07 


54.9 


51.0 


45.2 


B+- 


► K*+K^ 


V 


0.001 


0.001 


0.005 


-87.0 


-68.2 


519 


5°^ 




IV 


0.001 


0.001 


0.005 


-87.0 


-68.2 


519 




p^K"" 

r 


IV 


0.34 


0.35 


0.52 


-34.9 


-34.9 


-27.6 


B+- 


^p^K+ 


IV 


0.47 


0.41 


0.47 


19.0 


33.8 


51.5 


B^^ 


p-K+ 


I 


0.41 


0.47 


0.58 


-23.8 


-21.5 


-17.8 


B+ - 


>p+K^ 


IV 


0.02 


0.005 


0.05 


-84.4 


-90.4 


907 


B+ - 


► K*+'n 

1 


IV 


2.96 


2.95 


2.96 


21.6 


23.5 


27.3 


B+- 


> K*+r]' 


III 


0.27 


0.34 


0.86 


-59.3 


-5.79 


260 


B^^ 


K*^ri 


IV 


2.84 


3.13 


3.78 


22.6 


23.0 


23.5 


B^^ 




V 


0.08 


0.27 


1.23 


-75.7 


188 


292 


B^^ 




IV 


13.1 


14.7 


18.1 


52.5 


51.7 


50.2 


B^^ 




IV 


4.09 


4.93 


6.90 


67.8 


63.3 


56.1 


B+ - 




IV 


7.15 


8.02 


9.96 


44.5 


44.4 


44.1 


B+ - 


► K*07r+ 


IV 


11.5 


14.0 


19.9 


55.8 


51.7 


45.8 


B+- 


► (f)K+ 


V 


33.4 


17.9 


1.31 


51.2 


55.9 


118 


B°^ 




V 


31.6 


16.9 


1.24 


51.2 


55.9 


118 


B^^ 




V 


5.07 


0.01 


17.4 


53.1 


489 


31.0 


B+ - 




V 


5.31 


0.22 


21.2 


50.3 


-10.2 


28.7 
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The pattern B{r]K) < B{t]K*) < B{r]'K) and B{ri'K*) < B{r]'K) is found by CLEO |2g. 
For considered thirty seven B PV decays, three general features are as follows: 



• The theoretical predictions in the SM and TC2 model as given in tables are 
all consistent with the new experimental measurements and upper limits. 

• For most decay modes, the differences induced by using whether BSW or LQQSR 
form factors in calculations are small, ~ 15%. 

• The new electroweak penguin play a more important role for B PV decays than 
they do for B PP decays. 

For five B ^ pn and two B p'^rj^'^ decay modes, the NP contributions are very 
small, < 6% for N^-^-l" = 2 — cxo as shown in Table ^, and can be neglected. For B — > p^rj 
decay, the NP enhancement can be as large as ~ 110% for N^^^ = 3. 

For B uj-K decays, the NP contributions are small, < 13% for N^^^ = 2 — oo. For 
B ujri^"> decays, the NP contributions can be large but show a strong N^f^ dependence. 
The agreement between the theoretical prediction and CLEO measurement for B{B 
ujTT^) remain unchanged in TC2 model. 

For four B ^ (f)7i, (pr]^ ^ and four B — > K*K decay modes, the NP contributions can 
be as large as a factor of 4, but strongly depend on N^-^-^ . For two B (j)K decays, the 
NP enhancements are about 30% and insensitive to the variation of N^^^. It is clear that 
the Belle data of i? ^ K^(f) [^] prefer a small effective number of colors, say ~ 
For four Class-IV B —>■ K*tt decays, the NP enhancements can be as large as 90%, and 
are insensitive to the variation of N^^^ . 

For Class-I B^ — > p~K^ decay, the NP correction is about —20% and insensitive to 
N^^^ . For B^ p^K^ decay, however, the NP correction can be large in size, a factor 
of 17 enhancement for N^^-^" = oo, but very sensitive to the variation of N^f^. For the 
remaining two B p^K decays, the NP enhancements are large in size and insensitive 
to the value of N^^^ . 



4.3.1 B K*ri^'^ decays 

Very recently, CLEO reported their first observation |^ oi B ^ K*ri decays: 



-6 



B{B+ K*+7]) = (26.4+1^ ± 3.3) x 10 
B{B°^K*%) = (13.8t|^±1.6) X 10-^ 



while the theoretical predictions in the SM and TC2 model are 
B{B+ K*+7]) = 
B{B+^K*''ri) = 



1.5 - 3.8) X 10"^ in SM, 
1.9 - 6.1) X 10"^ in TC2, 



(1.5-4.5) X 10-6 in SM, 
(2.3 - 7.2) X 10-6 in TC2, 



(82) 
(83) 



(84) 
(85) 



where the uncertainties induced by using the BSW or LQQSE form factors, and setting 



;/2 ± 2GeV\ r] = 0.34 ± 0.08, N^^^ = 2 - oo, and = 200 ± 100 GeV, have 
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been taken into account. Although the central values of the theoretical predictions for 
B{B K*ri) decays are much smaller than the central values of the data, the theoretical 
predictions are still consistent with the data since the experimental errors are still rather 
large. Further refinement of the data will show whether there is a real difference between 
the data and theoretical predictions. The new physics enhancements to the decay rates 
are significant ( ~ 60% ) in size, insensitive to variation of N^-^-^ and hence helpful to 
improve the agreement between the theoretical predictions and the data, as illustrated in 



Figs.(HJ9) where the upper dots band shows the CLEO data |19, EO 



Fig.(H) and Fig.(|) show the mass and iV^'^'^-dependence of the decay rates B{B^ — >• 
K*^!]) and B{B'^ K*°r]), respectively. For Fig.||a and Fig.||a, we set N^'l^'f = 3. For 
Fig.||b and Fig.pD, we set m^- = 200 GeV. In these two figures, the dot-dashed line refers 
to the SM prediction, while the short-dashed ( solid ) curve corresponds to the theoretical 
prediction with the inclusion of NP effects from new gluonic ( both gluonic and electroweak 
) penguins. It is clear that the electroweak penguin play an important role for these two 
decay modes. 

For other two B K*!]' decays, the new physics enhancement can be significant in 
size, from —70% to ~ 200%, but strongly depend on the variation of N^^^ , as shown in 
Table |^. The theoretical predictions for these two decay modes are still far below the 
current CLEO upper limits. 



5 CP- violating asymmetries in B ^ PP, PV decays 

As is well-known, there are three possible manifestations of CP violation in B system [0, ^ 



55| , p6[ : the direct CP violation or CP violation in decay, the indirect CP violation or CP 
violation in mixing due to the interference between mixing amplitudes, and finally the CP 
violation in interference between decays with and without mixing. For the measurements 



of CP violation in B system, great progress has been achieved recently[^, 57 



In ref.[^, Ali et al. estimated the CP- violating asymmetries in charmless hadronic 
decays B — *■ PP, PV, VV, based on the effective Hamiltonian with generalized factor- 
ization. In another paper Chen et al. calculated the CP- violating asymmetries in 



charmless hadronic decays of Bg meson. We here will follow the same procedure of [26 1 to 
estimate the new physics effects on the CP- violating asymmetries of -B ^ PP, PV decays 
in the TC2 model. 

In TC2 model, no new weak phase has been introduced through the interactions 
involving new particles and hence the mechanism of CP violation in TC2 model is the same 
as in the SM. But the CP- violating asymmetries Acp may be changed by the inclusion 
of new physics contributions through the interference between the ordinary tree/penguin 
amplitudes in the SM and the new strong and electroweak penguin amplitudes in TC2 
model. The real and imaginary part of effective Wilson coefficients C^^^ and effective 
number will be modified by new physics effects and hence the pattern of Acp for 
two-body charmless hadronic B decays will be changed accordingly. 
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5.1 Formalism 



For charged B decays the direct CP violation is defined as 

in terms of partial decay widths. 

For neutral B^{B^) decays, the situation becomes complicated because of B^ — B^ 
mixing, and hence the time dependent CP asymmetry for the decays of states that were 
tagged as pure B^ or 5° at production is defined as 

Acp{t) - ^^^^^^^ ^ ^ ^^-^^^^ ^ . (87) 

According to the characteristics of the final states /, neutral B decays can be classified 



into four cases as described in For case-1, / or / is not a common final state of 

B^ and B^, and the CP- violating asymmetry is independent of time. We use Eg. (p6D to 
calculate the CP-violating asjTiimetries for CP-class-1 decays: the charged B and case-1 
neutral B decays. 

For CP-class-2 (class-3) B decays where B^ ^ {f = f) with = ±f ( /^^ ^ ±f 



) ||26[ , the time- dependent and time-integrated CP asymmetries are of the form 

•^cpit) = cte' cos(Am t) + ae+e' sin(Am t), (88) 
1 X 

^CP = — ntte' + — o(^e+e', (89) 

1 + X'^ 1 + 

where Am = rriH — rriL is the mass difference between mass eigenstates \Bfj > and \B^ >, 
X = Am/F 0.73 for the case of B^ — B^ mixing [Q, and 

1 - IAcpI' -2Im{\cp) 

> V,tVa <f\H,„\B«> 
- VW^</|ft„|B»>- 

For the formulae being used to calculate Acp for the more complicated CP-class-4 B 
decays, one can see Eqs. (36-40) of ref.p6|. We also define the ratio 

ATC2 _ ASM 
•^CP 

to measure the new physics effects on the SM predictions of Acp of B meson decays. 

As an example, we present here the explicit calculation for the Class-III-1 decay B^ 
7r^7r°. The decay amplitude M{B~ — »• 7r~7r'') has been given in Eq.(|59|) where all should 
be taken for transitions b —>■ d. For the charged conjugate amplitude we have 

M(E+^7r+7r°) = -t^UF,^-nml){ml-ml){V:,V^,ia, + a,) 

-VtlVtd X - (-ay + ag + aio + 08-^2) | , (93) 
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where the ratio R2 has been given in Eq.(^), and all are taken for transitions b 
The CP asymmetry for this decay mode is then defined as 



Acp{B' 



|M(5^ 



\M{B- 



\M{B- 



\M{B- 



TT TT 



0\|2 



71 TT 



0\\2 



d. 



(94) 



5.2 Numerical results 



In tables p|-p!4|, we present numerical results of Acp in 5 ^ PP and B PV decays in 
the SM and TC2 model. We show the numerical results for the case of using BSW form 
factors only since the form factor dependence is weak. In second column of Tables the 
roman number and arable number denotes the classification of the decays B PP, PV 
using iV^'^'^-dependence and the CP-class for each decay mode as defined in 



13 H, 

respectively. The first and second error of the theoretical predictions correspond to the 
uncertainties induced by setting Sk"^ = ±2GeV^ and Srj = ±0.08, respectively. 

The SM predictions for the CP-violating asymmetries of fifty seven B meson decay 
modes investigated here as given in tables are well consistent with those given in 
26|| . For details of the parametric dependence of Acp in the SM, one can see [|^. We 
here focus on the new physics effects on Acp of B meson decays. 

Very recently, CLEO reported their first measurements of CP-violating asymmetries 

for five decay modes[^, B^ K^tt^, Kgir^, uir^ and B^ K^tt^. They conclude 
that the measured asymmetries are consistent with zero in all five decay modes studied 
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Using the same input parameters as in Table |^, we find the theoretical predictions in 
TC2 model for those five decay modes 



AcpiB - 




(-3.4+J:^±0.8 


AcpiB ^ 




(-5.0t?;5± 1.1 


AcpiB - 




(-1.1 ±o.il[!:? 


Acp{B - 


.K^r,') = 


(-2.8lJ;^±0.7 


Acp{B 




(8.9 ± 0.11°:^ 



-0.7 +0.4 
-0.4 -0.3 
-0.1 +0.3 
"-0.2 -0.2 



1-2 



X 10" 

X 10-2 , 
±0.1 ±0.1) X 10' 



-0.8 
-0.5 
-1.7 
^11.1 



±0.1; 
±0.1) 



X 10 

X 10" 



-2 



(95) 
(96) 
(97) 
(98) 
(99) 

where the central values correspond to setting fc^ = ml/2, rj = 0.34 and N^-^^^ = 3, 
while the first to fourth error is induced by considering the uncertainty Sk"^ = ±2GeV'^, 
6r] = ±0.8, 2 < N;:^f < 00 and dm^ = ±100 GeV, respectively. For first four B ^ PP 
decay modes, the uncertainty induced by varying A'^^-^-^ is smaller or in comparable size 
with other three uncertainties. For B — >■ tuvr^ decay mode, however, the uncertainty 
induced by varying N^-^-^ dominate over other three uncertainties. 

The CLEO measurements, the 90%CL region and the theoretical predictions in the 



SM and TC2 model are as given in Table |T^. The theoretical predictions are taken from 
tables 9-14 and eqs.(^^-|^). One also should note that the sign convention as being used 
here in Eg. (|86|j87|) to define Acp is opposite to that used in [p^, we therefore changed 



the sign of the theoretical predictions of Acp in Table |I5] in order to be consistent with 
those reported results by CLEO. 

It is easy to see that the CP- violating asymmetries of all five decay modes studied are 
small in size in both the SM and TC2 model, and well consistent with the CLEO data. 
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Table 9: CP-violating asymmetries Acp in B 
p = 0.12 and N^^^ = 2, 3, oo for P = 
second error of the ratios corresponds to 6k 



PP decays (in percent) in the SM using 



J2 ± 2Ge1/2 and t] = 0.34 ± 0.8. The first and 
^ — ±2GeV'^ and drj = ±0.08, respectively. 



Channel 



Class 



oo 



5°- 
&^ 

&^ 

B^ 

B^ 

&^ 



TT TT 



rjrj 
rji]' 
ri'i]' 

> K^r] 



I- 2 

II- 2 

III- l 

II-2 
II-2 

II- 2 

III- l 

III- l 

V-2 

V-2 

IV- 1 

IV-1 
IV-1 

IV-2 
IV-1 
IV-1 

IV-2 

IV-2 
IV-1 

IV-2 



23.7 



+0.3+16.2 
-1.3-15.9 



-54.9: 



h9. 5+1.0 
^-3.9-1.2 
n 1 +0.01+0.02 
"J-J^ -0.02-0.01 

57^5+1.5+2.7 



-2.5-7.2 



60.811°: 



4.2 
0.7 



12.1^^ 



-0.3 
-0.2 
1 9 (^+2.9+0.9 
-'-^•'-'-6.3-1.2 

90 c:+2.8+5.2 
^0-3-1.5-5.8 

p;q cr+0.3+8.4 
'J'J-3-0.1-10.3 
f- f^+2.9+1.2 
"3-0-1.6-1.2 

1 +3.2+1.4 
~0--L-1.7-1.3 

-1.3 ±0.1 ±0.3 



34.4t[J:| 



-4-411:?: 



5.0 
6.4 

4 Q+1. 7+0.8 
^-U-3. 3-0.9 
9+1.0 
1.1 



34.7; 



-0.4H 
-0.6- 



5.0 
6.4 



-4.7 
-5.7 

1 n c;-t-5. 1+1.9 

-'-"-'-3_2.6-2.2 

1 o cr+S. 0+2.3 
-'-'-'-3-2.6-2.7 



90 /I +0.3+16.5 
^3-^-1.2-15.9 

1 c 0+6.1+2.6 
-'-3-3_3 4_3.o 

0.1 



hO.Ol- 
-0.02- 



0.02 
■0.01 



1 Q Q+4.7+2.3 
-'-3.o_2.6-2.8 

on n+5-9+3-5 
^U.U_3,3_4,o 



36.2tl| 



6.0 
7.2 

1/1 Q+2. 7+0.3 
-'-^-3_5.6-1.2 
1 p: cr+3.4+0.3 
±<J.O_7 2-0.9 

1^ q+5. 1+2.4 
-'-^-3_2.8-2.9 

94 Q+7.2+4.2 
/4.y_4 2_5.o 
f- Q+2.5+1.2 
~3.U_i 3_i 

9+3.2+1.4 
~0--^-1.8-1.3 

-1.2 ±0.1 ±0.3 



31.2 ±0.0: 



2-91^: 



5.9 



-0.7 
0.6 



-S.Qtli ± 0.8 
30.9 ± O.Otti 
31.2 ±0.0^^9 



10.4: 
13.4: 



-5.0+1.8 
-2.5-2.2 



90 n+0-3+16.6 
ZO.U_i 2-15.9 

48 n+1-0+7-8 

40-U-2.9-12.7 

0. 

f-q 1 +6.9+0.7 
~3-3.i_2.9_i, 5 

p;9 ^+6.9+0.4 
— OZ.O_2.9_i,5 

A J q+6.5+0.7 
-3-1.9-1.4 
1 o 1+2.7+2.6 
-Lo- -"--5.0-3.3 
99 /(+4.2+1.4 
^^-^-8.0-2.7 

Q Q+8.2+1.7 

-y.y_4.5_i.3 



-16.8 



+13+2.6 
-7.7-1.9 
Q Q+1-7+0.9 
~3.0_i 0-0.9 

-GAtli ± 1.4 
-1.2 ±0.1 ±0.3 

9P: ^+0.9+4.1 
-^3.0_o.6-5.0 

i.ol?:^±o.2 

-2.5li|:I ± 0.5 

90 7+1.2+3.9 
'^3- '-0.7-4.6 



qq 9+0.2 
33.^-0.3 

10.2^ 



-4.9H 
-2.7- 



2.2 
2.7 



5.0 
6.2 
5.0+1.8 
2.6-2.1 

1 q 1+4.9+2.3 
-'-3--L-2.6-2.6 



For all five decay modes, the new physics corrections are also small: which will change the 
SM predictions by about 20%. Fig.|10| and Fig.IlT] show the mass and N^-l"-^ dependence 



for B^ K^rj' and B^ ujir^ in the SM (the dots lines or curves) and TC2 model (the 
solid curves) Fl 



From the theoretical predictions and the CLEO measurements as given in tables ^[15 
the following general features about the CP-violating asymmetry of charmless hadronic 
B meson decays under study in this paper can be understood: 

• The CP-violating asymmetries depend weakly on whether we use the BSW or 
LQQSR form factors. The inclusion of NP contributions does not change this fea- 
ture. 



In these two figures we use the same sign convention as CLEO Collaboration l22| 



30 



Table 10: Same as in Table ^but in the TC2 model and assuming m^^ = 200 GeV. 



Channel 



B* 



TT ' TT 



rjr] 
rji]' 
rj'ri' 



Class 



I- 2 

II- 2 

III- 1 

11-2 
11-2 

II- 2 

III- 1 

III- 1 

V-2 

V-2 

IV- 1 

lV-1 
lV-1 

lV-2 
lV-1 
lV-1 

lV-2 

lV-2 
lV-1 

lV-2 



TC2 



27.3 

-55.5 
0.07 

51.9 

60.4 

55.3 
12.0 
13.1 

23.8 

46.2 
-3.8 

-4.8 
-1.2 

34.3 
4.11 
-3.34 

35.2 

30.1 
8.8 

11.5 



26.9 

14.9 
0.05 

10.7 

15.4 

27.5 
13.7 
15.9 

11.8 

21.3 
-3.4 

-5.0 
-1.1 

31.3 
2.7 
-2.8 

31.3 

31.2 
8.6 

11.3 



oo 



26.3 

49.8 
0.0 

-50.9 

-53.7 

-50.7 
15.7 
21.5 

-7.9 

-14.8 
-2.7 

-5.2 
-1.0 

26.1 
0.9 
-2.0 

24.7 

32.9 
8.3 

10.9 



CP\/0 



15.3 

1.0 
0.4 

-9.9 

-0.7 

23.3 
-0.7 
4.3 

-16.3 

-13.7 
-31.8 

-20.5 
-15.0 

-0.2 
3.5 
-24.6 

1.2 

1.1 
-16.4 

-15.2 



14.9 

-2.9 
0.4 

-22.4 

-23.2 

-24.0 
-4.0 
2.4 

-17.1 

-14.2 
-30.8 

-20.0 
-15.6 

0.2 
-3.6 
-22.4 

1.5 

0.04 
-17.0 

-15.8 



oo 



14.3 
3.7 

-4.2 

2.0 

7.2 
-13.5 
-4.1 

-19.8 

-12.3 
-28.9 

-19.1 
-16.4 

2.1 
-16.7 
-18.8 

4.1 

-1.0 
-18.0 

-16.7 



• The dependence of Acp is weak: 6Acp is about ±15% as one varies in the 
range lOOGeV < < 300GeV. 

• For twenty B PP decays, the new physics corrections to Acp are generally small 
or moderate in magnitude. The largest correction is about —30% for the decay 

K^Ti^, and about ±20% for decay modes B^ tt+tc- ,7:^7], K+7c~ , K^K^ 
and B^ K^tt^ , K^K^. For four class-11 B rjrj^'^rj'ri' and vr'^vr^ decays, they 
have large CP violating asymmetries ( around ±50%), but unfortunately also have 
strong A^^-^-^-dependence in both the SM and TC2 model. 

• For B —>■ PV decays, however, the NP corrections to Acp can be rather large 
for many decay modes, as illustrated in Table For class-1-4 decay B^/B^ 
p^7r^/p^7r+, the new physics correction is (60 ~ 100)% for N^-^-l" = 2 — oo. For 
decay B^ — > K*'^K^ the correction can even reaches a factor of 20 for A/"^-^-^ = 2 
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Table 11: CP-violating asymmetries Acp in B ^ PV decays (in percent) with b — > 
d transition in the SM using p = 0.12 and N^^'l' = 2,3, oo for k"^ = ml/2 ± 2 and 
rj — 0.34 ± 0.8. The first and second error of the ratios corresponds to Sk"^ — ±2 and 
577 = ±0.08, respectively. 



Channel 



Class 



-J-^-O.T 



1^ — 

18 

r 0+0.7+10.5 
^■°-1.8-8.8 

nc 1 +4.4+5.7 
<J'J-J--l.l-4.8 
A 1+2.9+0.8 
2-1.0 
9 7+0.7+0.4 
^- ' -1.5-0.3 

AQ. /I+10.7+2.3 
-4y.4_98_o.5 

o 0+0.5+2.7 
o-°+5.0-5.7 
A n+1-0+0.6 
3-0.5 
Q Q + 1. 1 + 0.9 
•5-^-2.4-0.6 

f-i 1+0.7+7.7 
'->^- ^-0.9-10.7 
1 r, 9+2.3+0.4 
-'-'-'■^-4.9-0.7 

r9 9+1.2+5.3 
'-'^■^-3.3-11.0 

-0.9+17.0 
-3.6-16.6 



q 9+1.2H 
"^•^-0.7- 



18 



00 



q 1+1.3+22.3 

"^•^-0.6-18 
cr 0+0.7+10.4 
^■0-1.8-8.9 

90 1 +0.4+16.9 
^"J-J^-l. 8-16.4 

I n 7+10.2+2.0 
"-■-U- < -4.6-3.1 

q (-+0.9+0.5 
<J-O-2.0-0.4 

0+2.3+1.8 
"J<J-0-4.6-4.5 

qi Q + 1.4+19.8 

'J^-y-6.5-17.4 
p- q+1.4+0.8 
^■^-3. 2-0.7 

f- Q + 1.6+1.1 

•J-^-S. 5-0.9 

qq n+0.1+11.6 
00. U_o.8-14.4 
9 -1 +1.6+0.4 
-0.6-0.6 

7 f^+0.+17.4 
' -0-0.3-14.6 

91 q+0.1+15.9 

^-'--'J-0.2-15.5 

II q+4.5+2.0 
-'--'--y -2.5-2.4 



(- 
5" ■ 



a;7r 

>■ CUTT 



± 



ujr] 
ujr]' 

077 
^ 077' 

B^B"^ 



K*^KyK*^Kl 



1-4 

I- 4 

II- 2 

III- l 
III-l 

II-2 

II- 2 

III- l 
III-l 

II- 2 

III- l 

II-2 

II-2 

V-2 
V-1 

V-2 

V-2 

IV- 1 

IV- 4 

V- 1 
IV-4 



32.3 



iD.U_4 j_3 2 
+5.6+2.4 
-3.0-2.6 

+5.6+2.7 
-3.1-3.2 

+5.6+2.7 
3.1-3.2 
+5.7+2.4 
-3.0-2.8 



12.7 

16.0 

16.0 
13.4 

13 3+5-8+2-5 



-7.9: 



-3.1-2.9 
5.0+1.2 

-17-0.8 
1 1 fi+3. 0+2.4 
-D-5.8-2.1 



r 0+0.7+10.4 

^■0-1.9-8.9 

214+8.6+3.5 

cr /I+3.9+1.0 
"J-^-LS-l.e 
q n+0.8+0.5 
<J-U-l. 7-0.4 

OA Q+7.7+4.2 

9(- Q+12. 5+3.4 
^"-^-39-2.0 
A cr+1. 1+0.7 
^-^-2.5-0.6 
4 g+1.2+1.0 



-1.7-0.6 

99 1+6.5+3.8 
^^--•--3.7-4.5 
o r+2.0+0.7 
o.O_44_o.9 

99 ^+5.4+3.9 
^^-^-3.0-4.6 



1 4+0-' 

-•--^-0.4- 



hO.2 
-0.3 

1 r,+0.8+0.2 
-'--'-'-0.3-0.1 

1 1+0.6+0.2 
-■--^-0.3-0.3 



1 4+o.6^ 

-■--^-0.3- 



hO.2 
-0.3 

19 ^+5.6+2.4 
-'-^-'-'-2.9-2.6 
1 q 9+5.7+2.4 
-'-<J-^-3.0-2.7 
A f^+6.3+0.5 
^-"-30.3-0.2 

-io-i;irk? 



9.1 



+4.7+1.6 

-2.4-1.9 



11.9 

11.9: 



h4.5+2.0 
-2.5-2.4 

1+4.5+2.0 
-2.5-2.4 
1 1 fi+5.3+2.1 
-L-L -0-2.8-2.1 
-r 9 q+5.5+2.2 
-'-^-<J-2.9-2.5 
7P; 0+0.6+8.6 
' 'J-O-12.2-12.3 
7 7+2.2+0.7 
~' -'-4.3-1.5 
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Table 12: Same as in Table O but with b —>■ s transition. 



Channel 



Class 



oo 



p-K+ 



K 



*± 



V 



K*^r] 

K* + TT- 



IV-1 
lV-1 

1-1 

lV-1 

lV-1 

III- 1 

IV- 1 

V- 1 

IV- 1 

lV-1 
lV-1 

IV- 1 

V- 1 

V-2 

V-2 
V-1 



1 c, 1+0.5+2.7 
-'-^•-'--l. 1-3.1 



+9.5+3.0 
-4.3-4.0 

+7.2+1.0 
2.9-0.7 
0.2+0.4 
-0.1-0.4 
7 0+4.1-1.2 
' •'J-2.1+1.5 
OO /1+16+3.6 
-/:y.4_2.4-1.9 



-17.9 

-11.7 
1.7 



5+0.6 
'-0.3 



±0.4 



-4.3t^i ± 1-0 

1 O+8.0+2.6 
10.0_4 4_2.0 

n 9+0.7+0.1 
U-^-1. 2-0.0 
1 1 Q+6.5+2.2 
-'^-'-•'J-3.6-1.8 

-1.5±0.ll[Ji 
-1.5±0.ll|]i 

31.1 iO.Olt^ 

oq I-+1.4+3.8 
^<J-3-0. 9+4.6 
1 1 p:+6.6+2.3 
-■^J^ •'-'-3.7-2.0 



32.1±0.0lii 

1 « Q+10. 1+1.9 

1 9 9+7.5+1.1 
-■-^•^-3.0-0.9 
9 7+0.6+0.7 
^- ' -0.3-0.6 
+4.1-1.3 
-2.2+1.5 



-7.3 



_c,4 4+22.0-0. 
'J^-^-0.8+2.4 



7 



-i.olj]:^ ± 0.2 

4 4+8.1+1.0 

^•^-1.1-1.0 

1 q Q+8.2+2.6 
10.J_4 5_2.0 

-1.7 ±0.0 ±0.4 
_1Q 4+6.0+2.1 
-'-U-'*-3.2-1.7 
1 I-+0.1+0.4 
-■-■^-0.1-0.3 

-1.6 ±0.1 ±0.4 
31.1 ±0.0^1^ 

oi 4+1.2+4.1 
<J-L.4-ll.4-5. 6 
1 7 O+10.3+2.7 
-i^»-4.1-3.6 



4g_0+1.4+1.2 
-9.7_2.7 

— 13 o+'^ °+-'^-^ 



0.6-3.6 
+5.1+1.8 
1.5 



-"-3.2-1.3 
+2.1+0.5 
'-7.9-0.6 
1+4.0+1.6 
--2.1-1.4 

-83.o;S:lii;| 



-2.3: 
-7.21 



0.1 
■0.2 



1 c; q+7.5+3.2 
10.0_]_3_3.4 



-14.0 



-8.2+1.6 
-4.6-1.9 



-4.3: 



}+l. 9+1.0 
-1.0-0.9 
o 7+4.9+1.9 
~0- ' -2.7-1.7 
-I 4+0.1+0.4 
-■-■^-0.0-0.3 

-2.5 ±0.1 ±0.6 

30.6 ± O.Oll^ 



24.2^1:1: 



4.0 
■4.8 



0.2l[5j±0.1 
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Table 13: CP- violating asymmetries Acp in i? — 
transitions in the TC2 model using p — 0.12, rj — 



■ PV decays (in percent) with b ^ d 
0.34, = ml/2, = 200 GeV and 



Channel 



Class 



TC2 



00 



SAcp[% 



00 



5°- 

B^ 

B^ - 



p~7l~^ 

p^T] 

p'^i 
> p^T] 

>p^r^' 



CUTT 







± 



> CUTT 

ujr] 
ujr]' 

H 



1-4 

I- 4 

II- 2 

III- 1 
III-l 

11-2 

II- 2 

III- l 
III-l 

II- 2 

III- l 

II-2 

II-2 

V-2 
V-1 

V-2 

V-2 

IV- 1 

IV- 4 

V- 1 

IV-4 



6.5 
7.5 

-36.7 
-4.3 
2.9 

-34.6 

-1.5 
4.2 
4.2 

48.8 
10.6 

56.2 

41.6 

16.8 
13.6 

16.8 

16.8 
10.6 

11.1 
-0.4 

-8.2 



6.1 
7.3 

20.8 
-5.8 
3.2 

16.8 

-16.6 
4.8 
4.8 

21.8 
8.9 

15.3 

66.0 

0.7 
0.4 

0.7 

0.7 
9.9 

10.5 
84.4 

-6.8 



5.4 
7.0 

24.7 
-11.0 
3.9 

58.4 

41.9 
6.3 
6.3 

41.0 
-2.2 

3.3 

23.1 

9.3 
6.9 

9.3 

9.3 
8.9 

9.7 
18.3 

-4.6 



104 
30.9 

1.7 
7.3 
6.6 

-29.9 

-117 
6.9 
7.7 

-4.5 
3.7 

7.8 

28.6 

5.4 
7.4 

5.4 

5.4 
-21.2 

-19.6 
-94.5 

-29.0 



88.1 
27.2 

-2.7 
7.0 
6.5 

-32.5 

-38.2 
6.7 
7.5 

-1.3 
5.0 

-31.5 

65.6 

-53.2 
-54.1 

-53.2 

-53.2 
-21.8 

-20.3 
-1951 

-33.4 



62.5 
21.1 

6.7 
3.1 
6.4 

-8.4 

20.2 
6.0 
7.1 

24.3 
8.0 

-55.9 

8.1 

-22.1 
-23.8 

-22.1 

-22.1 
-22.6 

-21.2 
-75.8 

-40.4 
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Table 14: Same as in Table O but with b ^ s transitions 



Channel 



B^ 
5± 



Class 



TC2 



oo 



6Acp\% 



oo 



IV- 1 
IV- 1 

I-l 
IV- 1 
IV- 1 

III- l 

IV- 1 

V- 1 

IV- 1 

IV- 1 
IV- 1 

IV- 1 

V- 1 

V-2 

V-2 
V-1 



19.2 
-8.9 



32.1 
-6.9 



1 o o 
— io.2 


And 

— 1 /.6 


2 8 


2 6 


—4 8 


— ^ 1 


-66.7 


-90.3 


-1.3 


-0.8 


-24.2 


-4.5 


-9.5 


-9.9 


0.2 


-1.6 


-6.6 


-6.2 


-1.2 


-1.2 


-1.3 


-1.4 


31.2 


31.2 


24.8 


31.5 


-8.5 


-19.9 



45.0 
-2.8 

-16.4 
-2.1 
-5.4 
-44.4 

0.3 

4.4 

-10.6 

-3.7 
-5.5 
-1.1 
-2.2 

30.8 

25.4 
0.06 



26.6 
-50.4 

55.2 
65.1 
-33.6 
127 

-27.2 

469 

-31.0 

0.5 
-41.6 
-18.9 
-13.2 

0.4 

5.8 
-25.9 



-0.09 
-63.3 

44.1 
-0.6 
-31.1 
65.9 

-20.9 

-203 

-28.5 

-8.0 
-39.8 
-19.5 
-13.2 

0.4 

0.4 
11.6 



-2.1 
-71.2 

26.0 
-9.5 
-25.2 
-46.5 

-52.1 

-71.2 

-24.0 

-14.3 
-36.4 
-20.3 
-13.5 

0.6 

5.1 
-66.9 



Table 15: CLEO measurements for Acp in 5 ^ Kti^Kti'^uott decays [22], and the 



corresponding theoretical predictions in the SM and TC2 model. 



Channel 


Aexp 
M-CP 


90% CL Region 




ATC2 
M-CP 


B^ K^n^ 

B^ UJTl^ 


-0.29 ±0.23 

-0.04 ±0.16 
±0.18 ±0.24 
-0.03 ±0.12 
-0.34 ±0.25 


[-0.67,0.09] 

[-0.30,0.22] 
[-0.22,0.56] 
[-0.17,0.23] 
[-0.75,0.07] 


[-0.001,0.079] 

[0.015,0.096] 
[0.007,0.017] 
[0.003,0.062] 
[-0.129,0.007] 


[0.009,0.058] 

[0.010,0.080] 
[0.006,0.015] 
[0.002,0.047] 
[-0.102,0.031] 
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due to strong interference between the contributions from the IV— penguin and new 
charged-scalar penguins. 

• For most class-I, III and IV decays, the iV^-^-^- dependence and /c^-dependence of 
6Acp are weak. For most class- V decays, however, the iV^-'^-^- dependence of SAcp 
is strong. 

• For most decay modes considered here, the new physics corrections on Acp in the 
TC2 model are still much smaller than existed theoretical uncertainties, and there- 
fore will be masked by the later. Low experimental statistics and large theoretical 
uncertainties together prevent us from testing the TC2 model through the studies 
of CP- violating asymmetries at present. 



According to relevant studies for these decay modes, we know that the FSI may 
provide a new strong phase and therefore enhance Acp to a level 20% — 40%, new physics 
with new large phases may also increase the Acp to a level 40% — 60%. Although there is 
still no evidence for direct CP violation in B system, the CLEO measurements ruled out 
large part of the parameter space for Acp- The key problem is that the measurements 
are currently statistics limited. 



6 Summary and discussions 

In this paper, we calculated the new physics contributions to the branching ratios and 
CP-violating asymmetries of two-body charmless hadronic B meson decays B PP, PV 
in the TC2 model by employing the NLO effective Hamiltonian with generalized factor- 
ization. We will present the calculation for the new physics effects on B{B — > VV) and 



AcpiB VV) in a forthcoming paper ||60 



We firstly evaluate analytically all strong and electroweak charged-scalar penguin di- 
agrams in the quark level processes b — >■ sV* with V = 'y,gluon, Z, extract out the 
corresponding Inami-Lim functions C^f*"^, -D^*"^, -Ej*"^ and Cj'""^ which describe the NP 
contributions to the Wilson coefficients Ci{M]v) {i = 3 — 10) and Cg{M]y), combine these 
new functions with their SM counterparts, run all Wilson coefficients from the high en- 
ergy scale fi ~ 0{Mw) down to the lower energy scale = 0{mh) by using the QCD 
renormalization equations, find the effective Wilson coefficients C^^^ , and finally calculate 
the branching ratios and CP-violating asymmetries after inclusion of NP contributions in 
the TC2 model. 

In section ^, we calculated the branching ratios for fifty seven B PP, PV decays 
in the SM and TC2 model, presented the numerical results in tables (^[]^) and displayed 
the rriji and A*"^-^ •'^-dependence of several interesting decay modes in Figs.(P-P). From these 
tables and figures, the following conclusions can be reached: 

• The theoretical predictions in the TC2 model for all fifty seven decay modes under 
study are well consistent with the experimental measurements and upper limits 
within one or two standard deviations. 
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• The theoretical uncertainties induced by varying fc^, rj and m^^ are moderate within 
the range of P = ml/2 ± 2GeV^, r] = 0.34 ± 0.08 and = 200 ± 100 GeV. The 
dependence on whether we use the BSW or LQSSR form factors are also weak. The 
Ai"^-^-^— dependence vary greatly for different decay modes. 

• For most B PP decay channels, the NP effects 63 are large in size and insensitive 
to the variations of the effective number of colors N^^^. For many B — * PV decays, 
however, 5B are sensitive to the variations of N^^^ . It seems that the B Kir and 
B —>■ Krj' decay channels are good places to test the TC2 model. 

• For most class-II, IV and V decay channels, such as -B ^ r/r/*-'^ B K-k,B Kt]', 
etc, the NP enhancements to the decay rates can be rather large, from 30% to 
100% w.r.t the SM predictions. So large enhancements will be measurable when 
enough B decay events accumulated at B factories in the forthcoming years. 

• For most decay modes, both new gluonic and electroweak penguins contribute ef- 
fectively. 

• For B —>■ Krj' decays, the new physics enhancements are significant, ~ 50%, and 
insensitive to the variations of /c^, r/, m^^ and N^^^ within considered parameter 
space. The theoretical predictions for B{B Ki]') become now consistent with the 
CLEO data due to the inclusion of new physics effects in the TC2 model. 

In section |, we calculated the CP- violating asymmetries Acp for B — > PP, PV decays 
in the SM and TC2 model, presented the numerical results in tables (|9|-|1^) and displayed 
the m^i- and A^'^-'^-^-dependence of Acp for decays B^ K^r]', ljtc^ in Figs.(|10|,|Tl]). In this 
paper, the possible effects of FSI on Acp are neglected. From those tables and figures, 
the following conclusions can be drawn: 

• Although there is no new weak phase introduced in TC2 model, the CP-violating 
asymmetries Acp can still be changed through the interference between the ordinary 
tree/penguin amplitudes in the SM and the new strong and electroweak penguin 
amplitudes in TC2 model. 

• The CP-violating asymmetries depend weakly on whether we use the BSW or 
LQQSR form factors in calculations in both the SM and TC2 model. 

• For most B —>■ PP decays, 6Acp are generally small or moderate in magnitude ( 
10% — 30% ), and insensitive to the variation of mfj- and N^-^^f . But the four class-II 

decay modes B^ tt^tt", r/*^'^//*^') have strong A"^-^ •'^-dependence in both the SM and 
TC2 model. 

• For B —>■ PV decays, however, 6Acp can be rather large for many decay modes. 
For decay 5°/^° — > p'^n / p , the new physics correction is (60 — 100)% for 
N^^^ = 2 — oo. For decay B^ K*^K^ the correction can even reaches a factor 
of 20 for N^^f = 2. For most class-I, III and IV decays, the A'^-'^-'-dependence and 
fc^-dependence of SAcp are weak. For most class- V decays, however, the N^^^- 
dependence of 5Acp is strong. 
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• For the measured five decay modes B — > Kit, Kr]' ,00%, the new physics effects is 
only about —20% w.r.t the SM predictions. The theoretical predictions for these 
five decay modes in the SM and TC2 model are well consistent with the CLEO 
measurements. 
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Appendix A: Input parameters 



In this appendix we present relevant input parameters. We use the same set of input pa- 
rameters for the quark masses, decay constants, Wolfenstein parameters and form factors 



as p|. 



Input parameters of electroweak and strong couphng constant, gauge boson masses, 
B meson masses, light meson masses, ■ ■ ■, are as follows (all masses in unit of GeV 

m 
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Mz 



1/128, as{Mz] 
91.187, Mw = 



0.118, sin^^ 



w 



^-5 



80.41, moo 

d 



0.135, nirj 



0.547, nirj' 



1.019, m^± = 0.494, m^^o 
1.64ps, r(5^) = 1.56ps, 



0.23, Gf = 1.16639 x 10 
= 5.279, m^± = 0.140, 
0.770, = 0.782, 
0.498, mK'± = 0.892, ttik^o = 0.896 



{GeVy 



mr,± = 
0.958, rrip 



(100) 



For the elements of CKM matrix, we use Wolfenstein parametrization, and fix the 
parameters A, \, p to their central values, A = 0.81, A = 0.2205, p = 0.12 and 
varying rj in the range of = 0.34 ± 0.08. 

We firstly treat the internal quark masses in the loops in connection with the func- 
tion G{xi, z) as constituent masses, 

mfc = 4.88^61/, nic = l.bGeV, m, = O.bGeV, rriu = ma = 0.2GeV. (101) 

Secondly, we will use the current quark masses for rrii {i = u, d, s, c, b) which appear 
through the equation of motion when working out the hadronic matrix elements. 



For p = 2.5GeV, one finds ||12 



mb = 4:.8SGeV, = 1.5GeV,m, = 0.122GeV, = 7.6MeV, m„ = 4.2MeF(102) 

For the mass of heavy top quark we also use mt = mlimt) = 168GeV. 

For the decay constants of light mesons, the following values will be used in the 
numerical calculations (in the units of MeV): 

A = 133, fK = 158, fK* = 214, /, = 210, ^ = 195, U = 233, 
/; = = 78, = = 68, = -0.9, f^, = -0.23, 

= -113, /,^, = 141, (103) 

where /^(/^ and /^^/j have been defined in the two- angle-mixing formalism with 6q = 
—9.1° and 6s = —22.2°^^ For more details about the mixings between i] and rj', 
one can see |61|, [1^ . 
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Appendix B: Form factors 



The form factors at the zero momentum transfer in the Baner, Stech and Wirbel 
(BSW) |T^ model have been collected in Table 2 of ref.|T2|. For the convenience of 
the reader we list them here: 

F(f-(0) = 0.33, F^^{0) = 0.38, Fo^''(0) = 0.145, Fq^^^O) = 0.135, 

<l2(0) = <i,2(0) = 0.28, <^*(0) = 0.32, AffiO) = 0.33, 

V^^iO) = V^^^iO) = 0.33, V^^'iO) = 0.37. (104) 

In the LQQSR approach, the form factors at zero momentum transfer being used 
in our numerical calculations are 



Fo^^(O) = 0.36, Fo^^(O) = 0.41, Fo^''(0) = 0.16, Fo^'''(0) = 0.145, 

{Ao, Ai, A2, V}{B ^ p) = {0.30, 0.27, 0.26, 0.35}, 

{Ao, A, A2, V}{B K*) = {0.39, 0.35, 0.34, 0.48}, 

{Ao, Ai, A2, V}{B ^uj) = {0.30, 0.27, 0.26, 0.35}. (105) 



The form factors Fo_i(/c^), Ao,i,2(^^) and V^k"^) were defined in ref. [|T5 

Fo{0) 



as 



Aoik') 

A2{k') 



1 - P/m2(0^ 

^o(O) 
1 - P/m2(0- 
^2(0) 



F,{e 

Ai{k 



l-k^/m?{l-y 

A(0) 
l-k^/m^{l+y 
V(0) 



(106) 



l-A;Vm2(l+)' ^ ' l-k^/m?{l-y 

The pole masses being used to evaluate the /c^.fiependence of form factors are, 

{m(0"),m(l'),m(l+),m(0+)} = {5.2789,5.3248,5.37,5.73} (107) 
for uh and dh currents. And 

{m(0"),m(l-),m(l+),m(0+)} = {5.3693,5.41,5.82,5.89} (108) 
for sh currents. 
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Figure 1: Typical self-energy and penguin diagrams for the quark level decays b {s, d)V* 
[V = 'j,Z'^,g), with (internal wave lines) and charged Pseudo-scalar exchanges (in- 
ternal dash lines) in the SM and TC2 model. The internal quarks are the upper type 
quark u, c and t. 
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Figure 2: Plots of branching ratios of S — > X+tt^ decay versus and l/N^^f in 
the SM and TC2 model For (a) and (b), we set N^^f = 2 and = 200GeV, 
respectively. The short-dashed line and solid curve show the branching ratio in the 
SM and TC2 model, respectively. The dots band corresponds to the CLEO data with 
2a errors: B{B K+tt^) = {ll-QtH) x 10"^ 
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decay mode. The dots band 
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Figure 5: Same as Fig.^ but for the case of 5 — K n decay mode. The dots band 
corresponds to the CLEO data with la error: B{B K°7r°) = (14.6^1^) x 10"^ 
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in the SM and TC2 model. For (a) and (b), we set N^^f = 3 and = 200GeV, 
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Figure 8: Plots of B{B^ —>■ K*^rj) versus m-jr and l/N^^^ in the SM and TC2 model. 
For (a) and (b), we set N^^^ — 3 and m^f = 200GeV, respectively. The dot-dashed 
line shows the SM prediction, while the short-dashed and solid curve refer to the 
ratios with the inclusion of contributions induced by new gluonic penguins and both 
new gluonic and electrowcak penguins, respectively. The upper band corresponds to 
the CLEO data with la error: B{B+ K*+ri) = (26.4;^°8^) x 10"^ 
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Figure 9: Same as Fig.(^, but for decay B{B^ K*^r]). The upper band corresponds 
to the CLEO data with la error: B{B° K*°r]) = (13.81^^) x 10"^. 
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Figure 10: Plots of CP-violating asymmetries Acp vs and l/N^^^ for decay 
{B^ — i> K^ri'). For (a) and (b) we set N^^^ = 3 and = 200 GeV, respectively. 
The MC.L. allowed region from CLEO is Acp = [-0.17, 0.23]. 
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Figure 11: Same as Fig.|T^ but for decay [B^ 
from CLEO is Acp = [-0.75,0.07]. 



UJ1T 



The 90%C.L. allowed region 



56 



